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Abstract: We study finite-size corrections to the magnon dispersion relation in three 
models which differ from string theory on AdSs x S 5 in their boundary conditions. Asymp- 
totically, this is accomplished by twisting the transfer matrix in a way which manifestly 
preserves integrability. In model I all world-sheet fields are periodic, whereas model II 
represents a particular orbifold of AdSs x S 5 and model III is a /3-deformed theory. For 
models I and II we construct the one-particle TBA equations and use them to determine 
the leading finite-size correction to the asymptotic Bethe equation. We also make some 
interesting observations concerning the quantization conditions for the momentum. For the 
same models we compute the leading and for model II the next-to-leading order finite-size 
corrections to the asymptotic magnon dispersion relation. Furthermore, we apply Liischer's 
formulae to compute the leading finite-size corrections in /3-deformed theory. In addition 
to reproducing known results, we provide new predictions for two-particle states from the 
sl(2) sector, to be confronted with explicit field-theoretic calculations in the dual gauge 
theory. Finally, we prove that the leading finite-size correction to the energy of an sl(2) 
magnon in orbifold theory is the same as the one for an su(2) magnon in /3-deformed theory 
for special values of (3. We also speculate that for these values of /3 our result for the next- 
to-leading order correction in the orbifold model might coincide with the corresponding 
correction to the energy of su(2) magnon in /3-deformed theory. 



*E-mail: G.E.Arutyunov@uu.nl, M.deLeeuw@uu.nl, S.J.vanTongeren@uu.nl 
t Correspondent fellow at Steklov Mathematical Institute, Moscow. 



Contents 



|]. Introduction 



[2]. Models: periodic, orbifold and /3-deformed 
|3|. Energy and Momentum 



I 







TBA equations for models I and II 

|4.1| Analytic properties of the asymptotic solution 

[4.2| The simplified TBA equations 



4.3| Exact Bethe equations 



Finite-size corrections to the dispersion: Model I 

l] LO wrapping correction to the energy and momentum 
.% Quantization of momentum and exact Bethe equation 



Finite-size corrections to the dispersion: Model II 
|6.1| LO wrapping correction to the energy and momentum 



6.2 Quantization of momentum and exact Bethe equation 
|6.3| Comparing energies of st(2) and su(2) states 
NLO wrapping correction to the energy 



Finite-size corrections in /3-deformed and physical orbifold theories 

|7.1| One-particle /3-deformed states 

|7.2| Two-particle /3-deformed states 

|7.3| Physical two-particle sl(2) orbifold states 

|7.4| Comparing /3-deformed and off-shell orbifold magnons 

Appendices 

T| Appendix A: Twisted transfer matrix 

Appendix B: Asymptotic critical points and TBA 
Appendix C: Wrapping corrections 
Model I 
Model II 



8.3.1 



8~3^ 



Appendix D: Details on the NLO computation 



8.4.1 The Harmonic number integral 



8.4.2 Some integrals 



12 

n 

17 

m 

24 
H 

m 

28 

m 



36 

36 
|8 
40 
40 

43 

n 

44 
46 
46 
48 
49 
49 





- 1 - 



1. Introduction 



In this paper we employ the mirror TBA in combination with adaptations of Liischer's 
perturbative approach to study single magnon configurations for various interesting models, 
which we discuss in detail below. This has been made possible by various important 
recent developments in the problem of finding the exact scaling dimensions of N = 4 SYM 
gauge-invariant composite operators with finite quantum numbers through the gauge-string 
correspondence [1]. 

First of all, Liischer's perturbative approach [2] has been adapted to the case of the non- 
Lorentz invariant AdSs x S 5 string sigma-model [3, 4], which led to a determination of the 
four- and five-loop anomalous dimensions of the Konishi operator [4, 5] and, subsequently, 
of all twist two operators [6-8] , see also [9] ; the corresponding results agree with first-hand 
field-theoretic computations [10-12]. 

Secondly, the Thermodynamic Bethe Ansatz (TBA) equations [13-17] for the AdSs x S 5 
mirror model 1 [19] have been constructed, based on the corresponding string hypothesis 
[20]. These equations should provide the tools needed to compute the energy of string 
states and, therefore, the scaling dimensions of the corresponding operators in the dual 
gauge theory, as has been nicely demonstrated on the example of the Konishi operator both 
at weak [21, 22] and intermediate [23, 24] values of the coupling constant. The computation 
of the five-loop anomalous dimension through the mirror TBA approach has been further 
extended [25] to all twist two operators and the results obtained agree with the ones 
based on the generalized Liischer's formulae [8]. Finally, we note that in general the TBA 
equations imply a set of functional equations known as the Y-system [26]. The Y-system 
for the case at hand was conjectured in [27], and its general solution was obtained in 
[28]. This Y-system has rather unusual analytic properties which are under investigation 
[13,17, 29,30]. 

In the present paper we continue to explore the mirror TBA approach. We concentrate 
on the TBA describing a one-particle excited state; a single magnon in the sl(2) sector 
with non-vanishing momentum p. This excitation does not exist in the spectrum of closed 
strings on AdSs x S 5 . The embedding fields describing closed strings must satisfy periodic 
boundary conditions which results in the level-matching condition; the total momentum 
of a state must be equal to zero. Of course, the concept of magnons remains of primary 
importance as it is used to define their scattering matrix; a fundamental building block of 
the whole integrability approach [31]. 

The spectral problem for a single magnon can be set up in different ways. One way 
is to start from the string sigma-model and require the embedding fields to be periodic. 
In this model, which we call model I, we can study more excitations upon relaxing the 
level-matching condition. Another possibility is to consider strings in other integrable 
backgrounds, where the single magnon is a physical excitation. Examples of this type are 
provided by strings on orbifolds of AdSs x S 5 and /3-deformed theories. 

lr The TBA approach in the AdS/CFT spectral problem was advocated in [18] where it was used to 
explain wrapping effects in gauge theory. 
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In general our interest in the single magnon spectral problem is three-fold. Firstly, 
we would like to understand the general relationship between energy and momentum (dis- 
persion) for integrable models in a finite volume and its dependence on the boundary 
conditions. Secondly, for theories where a single magnon is a physical excitation, it is in- 
teresting to compute finite-size corrections to the magnon energy, providing new data to be 
compared to the corresponding dual gauge theories. Thirdly, it is important to understand 
how to extend the TBA approach to a larger class of physical theories. 

In this paper we will analyze a one-particle excitation in three different models. The 
first is the already mentioned model I. Model II is defined by twisting the boundary con- 
ditions for sphere bosons by the total momentum of a state. As we will show, this model 
can be thought of as a particular orbifold of AdSs x S 5 . Finally, model III corresponds to 
/3-deformed theory. All models share the same world-sheet scattering matrix, which is also 
the one for strings on AdSs x S 5 [32-38]. As a consequence, the modified boundary condi- 
tions for models II and III in the large J limit 2 or at weak coupling can be implemented by 
means of the general theory of integrable twists [39]. Below we will discuss the modified 
boundary conditions starting from the underlying world-sheet theory. 

For models I and II we construct the corresponding TBA equations. Following the 
same strategy as [17], we assume that the TBA equations for excited states only differ 
from the ground state ones by the integration contours in the convolution terms. By 
requiring these equations to be compatible with the large J asymptotic solution given in 
terms of generalized Liischer formulae we then fix the shape of the integration contours. 
Simultaneously, we also determine how the TBA length parameter Ltba is related to J. 
Interestingly enough, for model I there is a noticeable difference between the level-matched 
and unmatched cases; Ltba = J + 2 in the first case and Ltba = J in the second. For 
model II we find an intermediate result for the length; Ltba = J + 1. 

Furthermore, by analytically continuing the TBA equation for the Yi-function to the 
kinematical region of string theory along the lines of [40] , we derive an explicit form of the 
exact Bethe equation Yi(u*) + 1 = for the rapidity variable u*. This equation should 
be viewed as a quantization condition for the particle momentum (rapidity) which replaces 
the asymptotic Bethe equation for the case of finite J. 

The large J asymptotic solution is constructed through Liischer formulae in terms of 
the transfer matrices in which the boundary conditions of interest are implemented through 
the corresponding twists. The Liischer formulae allow us to compute the leading order 
finite-size corrections to the energy Slo and momentum L\o . As was found for relativistic 
integrable models starting from an underlying lattice formulation, see example [41, 42], the 
energy E and momentum P are given by equations ( |3.1|) and ( [3.6D , respectively. In the case 
of non periodic boundary conditions, i.e. for models II and III, based on experience with 
the dual gauge theories, it is plausible to assume that the expression for the asymptotic 
energy is not modified. The situation regarding P is less clear however, and we just adopt 
formula ([O]) as a working assumption. We then show that for model I the leading finite- 

2 Throughout the paper J stands for the angular momentum of string rotating around the equator of the 
five-sphere. 
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size correction to the asymptotic Bethe equation coincides with Plo, while for model II 
this is not the case. A possible reconciliation of the two is that the introduction of a twist 
modifies the definition of the momentum in such a way that it becomes compatible with 
the exact Bethe equation. 

For model II, by employing some of the techniques of [5], we have computed the next- 
to-leading order correction to the energy -Enlo of an sl(2) magnon for J = 2,3,4 and 
6 3 . The results are given in terms of ("-functions and products thereof; all polygamma- 
functions arising in various steps of the computation are canceled in the final expressions. 
In these energy corrections, the term of maximal transcendentality has degree 2 J + 3. We 
also present results for the energy correction to some two-particle states from the sl(2) 
sector for a generic Z^-orbifold. As opposed to the single magnon, these states satisfy the 
level-matching condition, i.e. they are physical states of the orbifold model. 

Concerning /3-deformed theory, the construction of the TBA equations remains an open 
problem. Moreover, there is no consensus on the question of integrability of the finite-size 
model 4 . Here we restrict ourselves to the study of the leading Liischer corrections to the 
energies of states in the su(2) and sl(2) sectors. The Liischer formulae assume factorization 
of the full world-sheet S-matrix into a product of two su(2|2)-invariant S-matrices, and in 
this aspect they still reflect the properties of the large J solution. Understanding double 
wrapping should provide crucial insight on the issue of integrability of the /3-deformed 
theory in finite volume. 

Recently, wrapping energy corrections were found in /3-deformed SYM [44,45], which 
sparked attempts to reproduce these from finite size effects in the dual string theory. Close 
to what we consider here, an approach of obtaining the leading finite-size corrections to the 
energies of su(2) states in /3-deformed theories via twisting has been implemented in [46]. 
However, the twist we are using is slightly different from the one in [46] and we justify it 
from the world-sheet point of view. At the same time, our transfer matrix is also different 
from the one in [46], because it contains extra momentum-dependent prefactors needed to 
satisfy the generalized unitarity condition [19]. As a net result, the Yg-functions in the 
su(2) sector agree with the ones found in [46] for any number of excitations, and therefore 
in particular the Liischer corrections obtained there for a single su(2) magnon coincide 
with ours. Additionally, in [47] the /3-deformed analog of the Konishi state was studied by 
means of twisting certain elements of the S-matrix. Our twist incorporates this result, but 
derives it from a different starting point. 

We then turn our attention to the sl(2) sector in /3-deformed theory with /3 ^ j, where 
n is an integer. We determine the leading order corrections to the energy of two-particle 
states with J = 2,3,5 and give explicit formulae that can be readily used to compute 

3 Here the results take a nicely presentable form. 

4 For instance, the equations which encode a class of finite-gap solutions of the 7-deformed string sigma 
model have been discussed in [43]. However, these equations have been derived under an assumption 
that certain non-linear constraints are satisfied (in the limit 7 — » these constraints turn into the global 
conserved charges that do not belong to the Cartan subalgebra of SO(6)). It is unclear for the moment 
how to extend the corresponding approach to encompass all finite-gap solutions. We thank Sergey Frolov 
for the discussion of this point. 
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energy corrections for any J desired, providing new predictions to be confronted with 
future calculations in the dual gauge theory. We point out that for rational values /3 = j 
the correction starts at an order of g 2 higher than for generic /3. 

We also make a detailed comparison of orbifold and /3-deformed theories and show 
that for (3 = j the leading order finite-size correction to the energy of an sl(2) magnon 
in the orbifold model is the same as the corresponding quantity for an su(2) magnon in 
/3-deformed theory, in full agreement with the conjecture by [48]. We point out that if 
this interesting relation between the two theories continues to hold beyond the leading 
order, our result for Enlo found for the orbifold model will simultaneously give the NLO 
correction to the energy of an su(2) magnon in /3-deformed theory for these special values 
of /3. 

The paper is organized as follows. In the next section we introduce the three models and 
discuss implementation of the corresponding twists at the level of transfer matrices. Section 
3 contains some generalities on energy and momentum in the TBA approach. Section 4 is 
devoted to the construction of the TBA equations for models I and II. In sections 5 and 
6 we compute the leading finite-size corrections to the energy and momentum in models 
I and II, respectively. Section 6 also contains some results on the next-to-leading order 
correction to the energy of the st(2) magnon in model II. Section 7 contains our results for 
/3-deformed theories. Some definitions and technical details are relegated to appendices A, 
C and D. In appendix B we briefly discuss the issue of critical points for the asymptotic 
solution of the TBA equations. 

We note here that sections 2 and 3 are of general interest, which combined with section 
6 contain the background and explicit results on the leading order energy corrections in the 
orbifold theory, while combined with section 7 they contain the background and explicit 
results for /3-deformed theory and two particle st(2) orbifold states. Sections 4, 5 and 6 
contain the results which are important for next-to-leading order (and higher) corrections 
and our observations regarding momentum quantization. 

2. Models: periodic, orbifold and /3-deformed 

An important peculiarity of the string sigma-model on AdSs x S 5 is that the fermions 
exhibit unusual periodicity conditions related to the total momentum p [49,50]. Indeed, 
in the light-cone gauge one of the fields, an angle cp which parametrizes the five-sphere, is 
unphysical and must be solved in terms of physical (transversal) fields. The equation of 
motion for <fi then implies 

0(2tt)-#O)=p. (2.1) 

On the other hand, in order to make the world-sheet fermions neutral under the isometry 
<p —> <p + constant, and simultaneously bring the Wess-Zumino term of the string action 
to a local form, we need to perform a field redefinition of the fermions ip — > e^ip. As a 
result, the new fermions satisfy the twisted boundary conditions 

^(2tt) = e* p ip(0) . (2.2) 
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Since we aim to describe closed strings, all the sphere bosons must be periodic which leads 
to the level-matching condition p = 2irm, where m is the winding number. As a result, the 
theory splits into sectors with periodic and anti-periodic boundary conditions for fermions, 
depending on whether m is even or odd. 

Moving on, we recall that in the large J limit the symmetry algebra of the string sigma- 
model coincides with two copies of the centrally extended su(2|2) superalgebra with the 
same central charge depending on p [34, 51]. The fact that the 16 physical fields (8 bosons 
and 8 fermions) transform under two copies of su(2|2) suggests to treat this representation 
as a tensor product 4x4, where 4 is a fundamental four-dimensional representation of 
su(2|2) with two bosons and two fermions. 

In order to treat a single magnon, we have to relax the level-matching condition. It is 
then easy to see that it is impossible to twist the boundary conditions for factorized world- 
sheet fields in such a way as to have all sphere and AdS bosons strictly periodic, while 
fermions satisfying the twisted boundary conditions (|2.2|) . In the following we therefore 
introduce three models which have the same field content as the string sigma-model on 
AdSs x S 5 but differ in their boundary conditions. Within these models we will be able to 
study the spectral problem for a single magnon. 

Model I — Periodic fields 

In this model one assumes that all fields - bosons and fermions - are periodic. The level- 
matching condition is not imposed allowing to consider a single magnon with non- vanishing 
momentum. Clearly, the periodic boundary conditions preserve all the global symmetries 
of the original light-cone string sigma-model. 

Model II Orbifold 

In this model one twists the boundary conditions for bosonic fields y a and t/q, where a = 1, 2 
and a = 1, 2: 



2/i(2tt)\ = / \ / yi (0)\ / Vi&*)\ = / e<* \ / yi (0)\ 

y 2 (2v)) ^ e-*** J \y2(Q) J ' ^ 2 (2tt) ) { J \ y 2 (0) J ' 



while fermions 9 a and 9 a kept untwisted and so the AdS fields as well. Here we should 
distinguish two possibilities for choosing the left and right twists: 



The twists above will lead to a modification of the boundary conditions for two of the 
four fields Y a h- Together with the two fields contained in the complex Z field which is 
necessarily unaffected by the twist, these fields parametrize the five-sphere. Furthermore 
the physical fermions 9 aa = y a 9 a and r]aa = 9 a ya inherit the twisted periodicity conditions 



ae = a r 



or 



on = —a r . 



(2.3) 



9 la {2n) 
2d (27r) 



e iae 9 la (0) , 
e- iai 9 2a (0) 



T/ io (27r)=e^7/ la (0), 



(2.4) 
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The two twists in equation lead to essentially equivalent theories, so in what follows 
we have chosen to explicitly treat the case ct£ = —a r . 

A particularly interesting situation arises when a twist is chosen to be 

«<? = — , ( 2 -5) 

where n and S are integers. In this case the twisted boundary conditions coincide with 
the ones for string theory on the orbifold AdSs x S 5 /T, where T is a cyclic subgroup Zg in 
su(2) C 5u(4), see [52,53]. The canonical element generating Zg is ui = e~ and various 
values of n, n = 0, . . . , S — 1 describe the corresponding twisted sectors. Since two of the 
three complex scalars, Z and Kj, do not undergo twisting, the orbifold model preserves 
an su(2) part of the original su(4) R-symmetry on the bosonic fields. However, since all 
fermions are twisted, it has no supersymmetry 5 . 

Yet another possibility, which we will be mostly interested in, is to chose a? = —a r = ^, 
where P is the total momentum of a state. Excitations in this model with these boundary 
conditions are embedded in the orbifold model. Indeed, at large J the momentum of a 
single magnon is subjected to the quantization condition p = n G Z, which is nothing 
else but the Bethe equation, cf. our discussion at the end of this section. Therefore, for 
fixed J we deal with a magnon in the twisted sector corresponding to uj n = e~ of the 
orbifold theory. As we will see, the momentum p receives finite-size corrections starting 
at order g 2J+A ; these effects are compatible with the quantization condition, and hence 
preserve the orbifold structure. For LO and NLO corrections to the energy of a single 
magnon the difference between p and the total momentum will not play any role. 

In both model I and II we will study states from the s[(2) sector. For states from 
this sector the level-matching condition implies that the total momentum of a state must 
be equal to zero. In particular therefore a single s[(2) magnon is not a physical state of 
the orbifold theory. In what follows we will therefore study one-magnon configurations by 
relaxing the level-matching condition. In section 7 we also compute the leading finite-size 
corrections for a few physical two-particle states of a generic orbifold. 

It is worth pointing out that the orbifold boundary conditions with the twist parameter 
chosen to be proportional to the total momentum introduced here, are similar but not 
equivalent to those of /3-deformed theories. Nevertheless, as will show this model exhibits 
an interesting relation to /3-deformed theories. 



Model III — /3-defor mat ions 

This is a model which describes an exactly marginal deformation of AdSs x S 5 superstring 
theory. It has been introduced in [54] and further studied in [43,55-58]. The model 
admits a natural generalization to a more general non-supersymmetric background, which 
is obtained as a three-parameter deformation of AdSs x S 5 . The latter can be regarded 

5 We would like to thank Matteo Beccaria for correcting a statement about the presence of supersymmetry 
in this model in the previous version of this paper. 
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as the same AdSs x S 5 string sigma-model but with twisted boundary conditions for three 
isometric angles fa of S 5 , namely, 



fa{2ir) - fa(0) = -2vre ijfc7j J fc , i = 1,2,3. (2.6) 

In these 7-deformed theories, the j{ are three deformation parameters and the J{ are con- 
served Noether charges corresponding to the shift isometries in the fa direction. Because of 
the presence of the Jj, the deformations depend on the state of interest. For the symmetric 
case 

7i = 72 = 73 = P (2.7) 

with j3 real, the deformation preserves 8 supercharges and it is dual to the J\f = 1 supercon- 
formal field theory which is an exactly marginal deformation of N = 4 super Yang-Mills 
theory [54]. 

To continue, let us recall that the transverse bosonic fields describing the five-sphere 
can be combined in a matrix Y 

Y =( Y n Y iA = ( ae ~ i4>2 be ~ 1 * 3 ) (2 8 ) 

\Y 2i Y 22 ) \ be*** -ae^J ' 1 ' ; 

where we exhibit explicitly the dependence on two isometric angles fa and ^3 [50]. The 
other angle (p\ parametrizes the complex field 

1 _ n 2 - h 2 

Z = — -e l<t>l 

l + a 2 + b 2 ' 

Thus, the three </>'s together with a and b provide a parametrization of the five-sphere. We 
see that equation (^) implies the following twisted boundary conditions for sphere bosons 



e 



i7T/3(J 2 -J 3 ) Q \ / e i7T/3(2Jl-J 2 -J 3 ) Q 



Y ^ = [ o e -*r/j(*-j.) ) [ e ^mH-H-i,) ) ■ ( 2 - 9 ) 

Obviously, these conditions are inherited from the following periodicity conditions for y a 
(left bosons) and y„ (right bosons) 



[ |/i(2tt) ) 




V 2/2(2tt) J 




( yi(2vr^ 


H 


1%(27t) 1 





e i7T/3(J 2 -J 3 ) \ / yi ( ) 

e -i7r/3(J 2 -J 3 ) I I y 2 (Q) 



i7r/3(2Ji-J 2 -J 3 ) \ /j/i(0) 

e _^ ( 2J 1 -J 2 -J 3 ) J I .(g) 



(2.10) 
(2.11) 



Also, we have 

Z(2vr) = Z(0)e 2 ^( j2 " j3 ) . (2.12) 

As was mentioned above, in the light-cone gauge the field <f> = fa is non-dynamical and 
satisfies equation (|2.1[) . Thus, p = 2irf3(J2 — J3) is the level-matching condition in /3- 
deformed theory. Finally, we mention that the fermions 9 a and 9^ have standard periodic 
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boundary conditions, while the physical fermions 9 a a = y a 9a and rj a a = 9 a yd will inherit 
the twisted boundary conditions from the bosons. 

In /3-deformed theory we will be interested in two kind of states which we call states 
from the 5u(2) and sl(2) sectors respectively. The states from su(2) have J3 = and 
are characterized by two charges J = J\ and M = J2; the level-matching condition is 
p = 2it(3M. For J large such a state can be represented as 

A+( m )...A+(p M )|0), 

V v ' 

M 

where is a (bosonic) creation operator from the corresponding Faddeev-Zamolodchikov 
algebra [59]. The gauge theory operators dual to the states from su(2) have schematically 
the form ti(& M Z J ), where a complex scalar is dual to Y^. A single su(2) magnon with 
p = 2-7T/3 is a physical excitation of /3-deformed theory. 

The states from sl(2) have a single charge J and a spin S, and in the large J limit can 
be represented as 

A+.{p 1 )...A+.{p s )\0), 

V v ■ ' 

s 

where 3 and 3 are fermionic indices. Since J2 = M = 0, the level matching condition 6 is 
p = 0. A single sl(2) magnon can be considered off-shell, but in this case it is not a phys- 
ical excitation of /3-deformed theory. The dual gauge theory operators are schematically 
tr(D s Z J ). 

Twisted transfer matrix 

In the large J limit the spectrum of the models discussed above is captured by the asymp- 
totic Bethe Ansatz equations (the Bethe-Yang equations) which are constructed taking 
into account the corresponding boundary conditions. In the framework of the Algebraic 
Bethe Ansatz, (modified) boundary conditions are accounted for by means of twisting the 
transfer matrix. Now consider M string theory particles characterized by the u*-plane ra- 
pidities ui, . . . , um or, equivalently, with momenta pi, ■ ■ ■ ,Pm- Consider also an auxiliary 
particle with rapidity v corresponding to a bound state representation ttq of su(2|2) with 
the bound state number Q. Scattering this auxiliary particle through iV particles gives rise 
to a monodromy matrix 

M 

T(v\u) = Y\_^ai{v,Ui) . 

i=l 

Here §> ai (v,Ui) is the S-matrix which describes scattering of the auxiliary particle with 
another one with rapidity U{. As a matrix acting on the auxiliary space, T(v\u) satisfies 
the fundamental commutation relations 

§12(^1, v^T\(v\\u)T2{v2\u) = T 2 (v2\u)Ti(vi\u)E>i2(vi, v 2 ) ■ 
6 For simplicity we restrict ourselves to the case of no winding. 
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Sector 


sl(2) 


su(2) 


Factor 


L 


R 


L 


R 


Periodic (I) 














Orbifold (II) 


p/2 


±p/2 


p/2 


±p/2 


/3-deformed (III) 





2nJ(3 


vrM/3 


(2J-M)ir/3 



Table 1: The twists for the different models in the s[(2) and su(2) sectors. The off-shell su(2) 
magnon with zero twist has not been investigated. 



We can also introduce a twisted transfer matrix 

T{v\u) = Tr [ttqG/) T(«|tZ)] , (2.13) 

where an element 5 is a twist and trace is taken over the auxiliary space. If g is such that 
[§12,5 ® g] = 0, then the fundamental commutation relations imply that T(v\u) commute 
for different values of v and therefore define a set of commuting charges. For the string 
S-matrix we are interested in g £ SU(2) x SU(2) and as such the most general twist will 
involve six arbitrary parameters. Here the first SU(2) corresponds to the AdS space, while 
the second one corresponds to the five-sphere. 

Postponing development of the general theory to the future 7 , we note that for all cases 
of present interest the twist is of the form g = 1 ® K, where 

• For model I a = 0, i.e. twist is absent for both left and right copies of su(2|2). 

• For model II we have two options 

at = a r = - or ag = —a r = - , (2-15) 

for the left and a£ and a r refer to the left and the right transfer matrices, respectively. 

• For /3-deformed theory the left and right transfer matrices corresponding to left and 
right copies of su(2|2) are twisted with 

at = vr/3(J 2 - J 3 ) , a r = 7r/3(2J! - J 2 - J 3 ) (2.16) 

for the left and the right transfer matrices, respectively. 

For the reader's convenience we have summarized the twists for the different models in the 
st(2) and su(2) sectors in Table [TJ. 



7 A generic twist corresponds to describing the most general orbifold and 7-deformed theories; it breaks 
global symmetries of the model and brings the Bethe roots at infinity to finite positions in the complex 



plane [57,60]. 



- 10 - 



The Q-particle bound state representations of su(2|2) have dimension 4Q and are of 
two types; symmetric and anti-symmetric. Hence we can define two transfer matrices, 
Ti jS and T 1 corresponding to symmetric and anti-symmetric bound state representations 
respectively. The relations between Ti jS and T a>1 are schematically depicted in Figure 
[l]. These transfer matrices are now used to construct the asymptotic Y-functions Yq in 
the 5l(2) and su(2) sectors 8 . The corresponding expressions are given by the generalized 
Liischer formula [4]; 

M 

Y^{v) = e- J ^T\v\u)T r (v\u)\{S^ ) {v,u, i ). (2.17) 

i=\ 

Here £q is the energy of a mirror Q-particle, S®^{v, u) denotes the S-matrix with the first 
and second argument in the mirror and string regions, and T^ r is the left, respectively 
right twisted su(2|2) transfer matrix. 

Below we discuss the asymptotic Bethe equations which in the present approach arise 
from the condition Y\*{uj) = —1, the latter being the equations for the string rapidities Uj. 
Here Yu is the asymptotic Y\(v) function analytically continued to the kinematic region 
of string theory. 

Bethe equations in the st(2) sector 

In appendix A we present the full eigenvalues of the twisted transfer matrix Tq^ evaluated 
on the s[(2) vacuum, as well as the corresponding auxiliary Bethe equations which can be 
easily obtained by generalizing the result of [61] found for the untwisted case. For states 
from the s[(2) sector we have T e = T r = Tq^\ where the eigenvalue Tqi does not involve 
auxiliary Bethe roots and is given by 

tf'n^i = i+ii ; vi — 




2ik 
9 



«i + f(Q-l) 



Q-l M 
m=± k=l i=l 



Definitions of various quantities entering the last formula can be found in appendix A; cos a 
is a twist of the bosonic eigenvalues. The subscript of the transfer matrix indicates the 
anti-symmetric bound state representation which an auxiliary particle transforms in upon 
continuation of its spectral parameter from the mirror to string theory kinematical region. 
The transfer matrix is normalized as ]P = 1, so that the asymptotic Bethe equations in 
the s[(2) sector have the form 

M 

l = e^ J HS sim ( Uj ,u t ). (2.19) 



8 The reader should not confuse the asymptotic Yq- and exact YQ-functions. Since we are mainly in- 
terested here in the asymptotic solution, we will often omit the superscript "o", for instance, in section 4 
where the TBA equations are discussed. 
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Here S^)( u jj u i) = Ssl(2)( u ji u i) 1S the usual sl(2) sector S-matrix. In particular, for a 
single excitation this gives 

1 - Jp j 



From here, note that in none of the models we consider the twist enters the sl(2) Bethe 
equations. This of course should be the case because in all our models, including the 
/3-deformed one, the sl(2) bosons have usual periodic boundary conditions. 

Bethe equations in the su(2) sector 

In the su(2) sector the basic object needed to build up the Yq functions through the 
generalized Liischer formulae is the transfer matrix Tq^\ the left and right copies of 
which are generically twisted differently. For the twisted one has 

M _ 



= (q + 1) n 

i=i 

M 

-Qe ia H 

i=i 



M _ 
X 



X -1 



v - X„- X, 



X: V X " ±J - 



x — X 




X^X- 



X X 



(2.20) 



x n 



x+xt 



X X- 



+ (q - 1) n 



M _ 

X — X, 



i=l 



x n 



X, X 




This time the subscript of the transfer matrix indicates the symmetric bound state represen- 
tation which an auxiliary particle transforms in upon continuation of its spectral parameter 
from the mirror to string theory kinematical region. The asymptotic Bethe equations now 
read 



M 



l = Y u (u k ) = e^ J e 



iPkJ p -i{a l +a r ) 



x k -x 



i=i 



„+ 



In particular, for a single excitation we have 9 



i(pJ-a e -a r ) 



(2.21) 



(2.22) 



3. Energy and Momentum 

In the TBA approach for the AdSs x S 5 mirror model, the energy of an M-particle state 
from the sl(2)-sector is given by 




Here the integration runs over a real rapidity line of the mirror theory and fH are momenta 
of the mirror Q-particles. Moreover, £(p) is the asymptotic energy of a string theory particle 
with momentum p, given by 

£(p) = ^1 + 402 sin 2 !- (3.2) 

9 Recall that S , s( ( 2 )(u,u) = — 1. 
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Figure 1: The relations between eigenvalues of the transfer matrix calculated for (anti-)symmetric 
representations on a bosonic (u>) and fcrmionic (9) vacuum [62]. BR stands for the Bazhanov- 
Reshetikhin formula [63], c.c. for complex conjugation or equivalently for switching the scattering 
matrix S for S^ 1 . 



The last term in the formula (3.1) can be understood as the finite-size correction to the 
asymptotic, i.e. large J, dispersion relation. 

As in the general case, for a single excitation the value of p, or equivalently the value 
of the rapidity u(p) 



u{p) = I C ot|yi + 4 5 2 sin 2 |, (3.3) 

is found from the exact Bethe equation Yi*(u) = — 1. As was pointed out in section 2, 
when J becomes large, the exact Bethe equations should turn into the asymptotic Bethe 
Ansatz equations, which for a single excitation are simply 



e 



ip J = 1 =>. p J = 2vrn , n = 0, 1, 



(3.4) 



Analogously, in various known integrable models the TBA approach leads to the fol- 
lowing expression for the total momentum 

— 1 00 r rlfQ 

p =E^-^E/ d ^Mi+^)) ! (3.5) 

i=l Q=l 

where here £® would be the energy of a mirror Q-particle; the integral term in the large 
formula can be viewed as a finite-size correction to the total momentum of string theory 
particles. For a single magnon one gets 

1 00 r df® 



where p is determined by equations (|3.4| ). It might appear slightly surprising that even for 
a single particle its momentum is modified by the finite-size correction. One might offer a 
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physical interpretation of P as being the combined momentum of a string particle and the 
mirror background. The correction term would then be thought of as the back-reaction 
of the mirror background to the particle; the former also gets into motion when a string 
theory particle is excited. 

On general grounds we could expect that in a finite-size theory it is the momentum 
P which must be quantized. Thus, we are led to two apparently different quantization 
conditions 

e iPJ = l ^ Y u {u) = -1 (3.7) 

which nevertheless should be compatible with each other. For model I we will find pertur- 
bative evidence that this is indeed the fact. However, in model II we will show that with 
the same definition of P only one of the two quantization conditions holds, namely the 
second one. 

4. TBA equations for models I and II 

In this section we start by discussing the analytic properties of the asymptotic solution and 
then show how these properties are taken into account in the excited state TBA equations. 
Finally, we summarize the complete set of excited state TBA equations for a single magnon 
in both the twisted and untwisted cases. These equations will be subsequently used to find 
NLO corrections to the energy and momentum through the TBA approach. 

4.1 Analytic properties of the asymptotic solution 

The analytic properties of the asymptotic Y-functions play an essential role in constructing 
the excited state TBA equations. This is because the integration contours for the excited 
state TBA equations are not known a priori. Demanding the TBA equations to be compat- 
ible with the asymptotic solution provides the necessary information to fix the integration 
contours 10 . 

In general we could expect that with the level-matching condition relaxed, the analytic 
properties of the asymptotic solution would become more intricate, and as it turns out this 
is indeed the case for the models of interest here. In addition we will see that the analytic 
properties of the corresponding solutions in models I and II are considerably different. Now 
as the TBA equations essentially involve three types of terms 

logy, log(l±y) and log(l-Y _1 ), (4.1) 

this means that as far as integration along a contour in the complex plane is concerned, 
we are interested in those points where 

Y(r) = 0,oo, ±1. (4.2) 

Of course, this approach is based on the assumption that the exact solution develops no new singularities 
in addition to those which are already present in the asymptotic solution. 
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When changing either the coupling constant g or the momentum p of a state these points r 
move in the complex plane, and as we will see, these points can either always lie on one side 
of the integration contour, or cross it dynamically when changing either g or p. Because 
the next sections of the paper will be directly concerned with the small g behaviour of the 
TBA equations, we relegate a discussion of so-called critical behaviour in g to appendix B, 
focusing here on the points r that are relevant already at weak coupling. 

Before moving on to discussing the particular zeros and poles in either model, we want 
to briefly comment on how zeroes and poles of the asymptotic Y-functions and points where 
they equal ±1 are related, and how to relate these to zeroes and poles of the transfer ma- 
trices; we will illustrate this by using the Y^i^-functions as a simple example. Analogous 
relationships hold for other Y-functions as well. 

In both model I and model II, because the denominator remains finite, zeroes of Y M \ VW 
are inherited from zeroes of the transfer matrices used to construct them 11 . The presently 
relevant Y-functions are constructed in terms of the transfer matrices as follows; 

l,M 1 l,M+2 v _ J-2,1 v _ -t2,3 J 2,l 



1 M\w — — ' 1 - ~ ~~ ' 1 - 



J-2,M+1 1-1,2 J- 1,2 J- 3,2 

T m,iTm+2,i Tm,iTm+2, 

Tm+i,2 'Im +11 T m+11 — Tm,iTm+2,i 



Y M\vw — — Tf, — — 777+ m 77\ • V 4 - 3 ) 



Here the ± superscript denotes the usual shift of the argument by ±i/g, and of course all 
transfer matrices are implicitly the sl(2) ones. From the last expression we immediately 
see that a zero r of Tm,i is also a zero of Yj^i vw , and it is clear that shifting this zero by 
±i/g we get the value of the rapidity where a lower Y-function becomes equal to minus 
one; Y^-i^O^) = -1, where r ± = r± i/g. 

While the construction of Ym\ w in terms of transfer matrices does not immediately 
allow for a simple interpretation like the one above, it should come as no surprise that 
similarly a zero of Y M \ W can give rise to points where Y M _i\ w is equal to —1. These 
relations for Y M \ W and Y M \ VW have been summarized in the diagram in Figure ^. Here we 
will not be concerned with poles, as they play no role in either model at small coupling. 

With these ideas in mind, we are left with a rather elaborate discussion of all zeroes 
and poles of the asymptotic Y-functions in both models. 

• Model I 

In the untwisted case, the Y^^-functions exhibit zeroes, while Y + and Y_ exhibit both 
zeroes and poles that have to be taken into account already at weak coupling. As men- 
tioned earlier, these zeroes are directly linked to zeroes of the transfer matrices; here each 
transfer matrix Tm has a single zero, which we will denote by tm—i- These zeroes have 
a natural ordering as tm > ^m+i and they give rise to two zeroes of Y M \ VW , tm-i and 
rjvf+i, and in addition to Y M \ vw (rj^) = —1. Also r\ is a zero of both Y + and Y_, as also 



The asymptotic Y-functions are constructed in terms of transfer matrices based on the underlying 
symmetry group of the model [64, 65] ; in the context of the string sigma model the corresponding asymptotic 
Y-functions were presented in [27]. 



- 15 - 




Figure 2: Zeroes of Ym r v \ w (+1) and their interrelations. Note that Tm here stands for the full 
set of zeroes of Tm+i- 



the transfer matrices in their denominators remain finite at this point. Finally, at u the 
numerators and denominators of Y± have poles. However in Y_ they cancel, whereas in Y + 
the numerator has a pole of higher order, giving rise to poles of Y + at . 

Furthermore Y + = 1 at the points ro and its conjugate fo- These points give rise to a 
string of poles of Y M \ W and Y M \ VW , which are located at tm and rjf, where ym = ro + iM/g. 
These point play no role at weak coupling however; they are discussed in more detail in 
appendix B. 

• Model II 

As mentioned above, twisting significantly modifies the analytic properties of the asymp- 
totic solution. First of all both Y M \ W and Y M \ VW now have zeroes. In the case of Y M \ VW , 
according to our general discussion, these zeroes originate from the ones of the transfer 
matrices; this time for M > 1 each transfer matrix has two zeroes, which are real or purely 
imaginary, and come with opposite signs. Explicitly, we have 

T M)1 (±r M -i) = 0, (4.4) 

where for M = 1 we only have T\ i(— ro) = 0. In fact, we were able to guess the exact 
analytic form of these zeroes. It looks very intriguing 

r M = V u2 - M{M + 2)£ 2 /g 2 , (4.5) 

where u is the Bethe root and £ is nothing but the asymptotic energy ( |3.1| ); note that 
ro = u. For p > 0, the last formula can be rewritten as 

r M = ^p ,/l-(M + l)2 s in 2 g. (4.6) 
g sin | V 2 



- 16 - 



In general we thus have Y M \ vw (±rM±i) = 0. For M = 1, while the story is a little different, 
the conclusion remains the same; Y\\ vw still has u as a zero although it is not a zero of T%, 
this is due to the pole of T2 at u~ . Of course, these zeroes give rise to points r M where 

^M|OTo(± r Af) = _ 1- 

Moving on, all Y^i^-functions turn out to have two zeroes which are always real. We 
will denote these zeroes by pm-1 and Pm+i) 

Y M \w{PM-l) = Y M \ w (p M +i) = 0. (4.7) 
Again we were able to guess their form, which fascinatingly enough is given by 

p M = -u M +i (4.8) 

where 

u M = f cot + 4(^)2 sin 2 (4.9) 

where one immediately recognizes that po = —u. Similar to Y M \ VW , 1 + Y M \ W has zeroes at 
the points pu ± i/g- 

Finally the Y± functions each have the same two zeroes 12 , Y+(±ri) = Y_(±ri) = 0. 
The points ±r^ f which gave rise to roots of 1 + Y M \ VW , here for M = split nicely into 
zeroes of 1 — Y- at —u^ and poles of Y+ at giving a very symmetric situation, as 
illustrated in Figure [|. In addition both Y+ and Y_ have a pole at p\ . 




Figure 3: Special points for Y± in the complex v plane. Note that r± can also be imaginary as 
opposed to real as is illustrated here. 

For the reader's convenience, an overview of the relevant roots and poles in both models 
is presented in Table ^. 

4.2 The simplified TBA equations 

In order to obtain the excited state TBA equations, following [17], we employ the con- 
tour deformation trick. In short we assume that the ground state and excited state TBA 

12 Not surprising since Y+ and Y_ share a common factor, which then of course is the source of the zeroes. 
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Y-function 


Zeros I 


Poles I 


Zeros II 


Poles II 


Cuts 


Y, 




?/± 
(X 




it 1 * 1 , pi 




Y, 1 


r ; ^0 










Y 


' 1 




— 1 — 


Pi 


n +9 


y_ - 1 






-«± 










r M j r m 


Pm±i 




±M,±(M + 2) 


Ym\w + 1 






Pm 






*M\vw 


TM±1 


r m , f m 


±r Af ±i 




±M,±(M + 2) 


Ym\vw + 1 


-± 
' M 











Table 2: Essential analytic structure of the asymptotic solution for a single magnon in both the 
twisted (right) and untwisted (left) case. A dash indicates absence of a zero or pole, while the empty 
spots would have information that is already in the table. The branch cuts are (— oo,2] U [2,oo), 
shifted up and down in the complex plane by the multiples of i/g indicated. N.B. —u^ should be 
read as — u ± ^ everywhere. 

equations differ only by the choice of the integration contours, and upon deforming the in- 
tegration contours of the excited state TBA equations to the ground state ones, we pick up 
contributions of extra singularities, leading to the appearance of new driving terms in the 
excited state TBA equations. This procedure is very similar in spirit to the Dorey-Tateo 
trick [40]. Below we discuss the choice of the integration contours in detail. 

The integration contours for all functions Yq, Ym\ w an d Ym\vw can be taken the same, 
hence we will speak of the integration contour for Yq, Y m \ w and Y M \ VW functions. This 
contour needs to be chosen such that it encompasses the pole of log 1 + Yq at the Bethe 
roots in the string region, which is naturally motivated by the fact that this gives a driving 
term for each Bethe root, see [17] for an explicit discussion. With this requirement it is 
only natural that the integration contour should run a little bit below the line — i/g after 
deforming it, and in fact we find that accepting or discarding contributions of zeroes and 
poles according to this choice of the contour, we obtain a set of TBA equations which is 
perfectly satisfied by the asymptotic solution. Returning the integration contour to the 
real mirror line, we pick up a contribution from all zeroes and poles that lie between the 
real line and the line —i/g, as well as the usual Bethe root contributions. The integration 
contour and the extra contributions that arise by returning the integration contour for 
excited states, C exc , to the real line, C can , are illustrated in the left panel of Figure § 

Analogously, for Y + and Y_ the integration contour needs to be picked essentially in 
the same manner, with of course the only exception that it should start at —2 and end at 
2. In particular, it is chosen in such a way that we pick up both the pole of Y + at u~ and 
the zero of 1 — Y_ at — u~ , as illustrated in the right panel of Figure |j. 

Now we would like to briefly come back to the zeroes of Ym\vw m m odel II. Looking 
at equation (|4.5| ), we see that they are either real or imaginary depending on a state (p) 
and the number M under consideration. Also, considered as a function of momentum, a 
given zero will go from real to imaginary values as the momentum is increased (0 < p < ir) . 
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Figure 4: The integration contour for all Y-functions except Y± (left) and the integration contour 
for Y± (right). 



When it is real or imaginary but lies between —i/g and the real axis, we need to take its 
contribution into account in the TBA equations. However, from equation (^1]) we see that 
as momentum is increased, a given zero will move out of the physical strip — l/g< Im < 1/g 
and eventually no longer needs to be taken into account. So the appearance of the zeroes 
rj\j in the TBA is state-dependent. Of course, the zeroes also move as g changes, but for 
now we are concerned only with the small g behaviour just discussed; we refer the reader 
to appendix B for further details on the behaviour in g. 

We would like to note that the kernels and S-matrices appearing below have been 
defined and are completely listed in [17]; for the sake of brevity we refrain from reproducing 
the list here and hence refer the reader there for their exact form. 

Taking the integration contours as indicated above together with the simplified TBA 
equations for the ground state derived in [13, 14], and returning the integration contours 
to the real line of the mirror theory, we obtain the excited state TBA equations for a single 
magnon. For both comparison purposes and brevity, we will discuss model I and model II 
side by side, splitting the TBA into the always present ground state and Y\ contributions, 
and the extra terms that come from zeroes and poles discussed above which differ between 
the twisted and untwisted case, see Table [|. Please note that in what follows for model II, 
for completeness we are always writing down the contributions from zeroes tm; of course 
they are not present if the corresponding zero falls outside the lower-half of the physical 
strip. We then find the following sets of integral equations. 

• M|w-strings; M > 1 , Y |i« = 

V M \w = logF M |t« - log(l + Y M _ 1 \ W )(1 + Y M+1 \ w )*s- 8 mi log _ j_ * s , 

1 y+ 

where for the two different models we have 

'o I 



Vm\w 



-logS'(p M _ 1 )5'(p M+1 ) II 



(4.10) 



The extra contributions in model II arise from the zeroes of 1 + Ym—x\ w and 1 + Ym+i\wi 
as indicated by the blue underlining; color coding has been used in all equations below for 
the readers convenience. Therefore, we see that while for model I these equations coincide 
with the ground state ones, for model II they do not. 

• M|fu>-strings; M > 1 , Y \ vw = 



Vm\vw= logY^i^^) 

+ log(l + Y M +i) *s- log(l + + Y M+1 \ VW ) -k s 

-S M i log \ — — * s , 
i - y + 



where for the two different models we have 



V, 



M\vw 



- S M i1ogS{u -v) - log5(r M _ 1 - v)S(r M+1 - v) 
- 5 M ilogS{±u~ - v) - logS'(±r Af _ 1 - v)S(±r^ +1 - v) 



II 



(4.11) 



Here and in what follows below we implicitly assume a sum over the ± sign in front of the 
zeroes. Given the shape of the integration contour, for M = 1 the first term in both models 
comes from the pole of 1+ at v = u~, and the zero of 1 — Y_ at v = —u~ in model II, 
whereas the second terms come from the zeroes of 1 + Y M _ 1 \ VW and 1 + Y M _^ VW , cf. Table 
|||. So we see that also at weak coupling, single particle states have additional contributions 
in the TBA equations for ?;u>-strings. Finally, we see that the pole of S(±u~ — v) cancels 



the zero of Y±\ vw {v) at v = ±u as it must, and that equation (|4.11| ) is compatible with the 
reality condition for Y- functions. 

• y-particles 



V± = log ^(v) + log S Uy (u, v) - log(l + Y Q ) * K Qy = , 
V T = log Y„Y + (v) + log i&- * S ( U) W ) _ 2 log 1 + Yl } vw * 



So 



1 + Y 1{1 



+ log (1 + Y Q ) * K Q + 2 log(l + Y Q ) * Kg 1 ★ s . 
where for the two different models we have 



-21og5(rf -v) I 
-21og5(±rf -v) - 21og5(pi ~ v ) 11 



(4.12) 
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Here the additional contributions come from the term log , , -J 1 ""' and we have used the 
following notation 

log * s(u, v) = r dt log f^j" s{t - V ) . 

#2 J~oc S 2 {U-t) 

Moreover, Si fy (uj,v) = Si y (z*j,v) is a shorthand notation for the S-matrix with the first 
and second arguments in the string and mirror regions, respectively. The same convention 
is used for other S-matrices and kernels. 

• Q-particles 

The TBA equations for Q-particles can be cast in particularly convenient form called 
the hybrid equations; their general form has been derived in [17]. It is here that we 
find another significant difference from level-matched cases studied previously, as well as 
between model I and II. We find that in order for the asymptotic Y-functions to satisfy 
the following integral equations for one-particle states, for model I the length parameter 
L TBA has to be equal to the angular momentum 13 J, whereas in model II we find that 
Ltba = >/ + 1. 

This merits a further discussion, as this behaviour is of course far from an accident. We 
have briefly investigated the hybrid equations for two and three level-unmatched magnons 
in the untwisted theory, and there we consistently find the relation to be Lj< BA = J. On the 
other hand, for level-matched states it was observed in [17] that L TBA = J + 2. Thus, in the 
untwisted case there is a clear difference between physical (level-matched) and unphysical 
(level-unmatched) states, rather than a difference between different numbers of magnons. 
Doing the twist removes some of the differences between the level-matched and unmatched 
cases, hence leaving us in some sense somewhere in between 14 with L TBA = J + 1. In the 
equations below we thus have 



J I level-unmatched 
J + l II 

J + 2 I level-matched 



The hybrid equations for Q-particles are now given by 

V Q = log Y Q (v) + L TBA £ Q - log (1 + Y QI ) * (k$$ -2s* K^' Q ) + log (u, v) 

- 2 log (1 + Y 1{vw ) * s * KyQ - 2 log (1 + Y Q _ 1{vw ) * s (4.13) 



2 log — is* Kl2 x - log r * K Q - log (1 - — ) (1 - — ) * p . v KyQ , 

1 — Y+ I — T7- I- 1 + 



13 In other words, the same as for the ground state [29]. 

14 We will see more examples of how the twist brings us 'closer' to the level-matched case later. 
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where for model I we have 



VQ 



2 log S*p. v . K y Q (r 1 ,v) - 21og5(r ( 



Q-i' 



v)+ 



(4.14) 



+ 2 log S -k p . v . K l v 2 x (u ,v)- log Sl$ x (u,v) 



while for model II we have 



Vq 



2log Si p . v K yQ {±r-,v) -21ogS(±r, 



Q-l' 



u)+ 



(4.15) 



+ 2 log 5 *,.„. ^(±u-, u) ± \ogS^ x (±u,v) + 



^og S Q (-u ,v) - logS yQ (-u ,v). 



In the above K V wx = and Yq\ vw = 0, which implies that the logS(±r^,z;) terms are 
also not present. The principal value prescriptions are required due to the pole of S(v) at 
v = —i/g, and are taken in accord with the integration contour C\ in the term log S-kK y Q 
running slightly above the real line, while for log S-kK^wx the contour C2 runs slightly above 
the real line in model I (like C±) and slightly below and above the negative and positive 
half real axis in model II respectively, as illustrated in Figure ||. We would like to note 
that the principal value prescription for the terms involving K y Q give extra contributions, 
but these precisely cancel against one another. All extra terms above follow from the same 
reasoning as for the other TBA equations; once more we have explicitly indicated their 
origin via color coding. 



Figure 5: The integration contours for the kernels log S * K v q and logS-kK^ x in the hybrid 
equations. 

4.3 Exact Bet he equations 

Analytically continuing the hybrid equations in the rapidity v to the string region needs to 
be carefully done as certain terms in the equations below have branch cuts when considered 



C, 




■2 



as a function of v. This continuation has been carried out in [17] for the Konishi state, and 
here we follow the same logic to obtain 



V Q , = vri(2n + 1) - iL TBA p - log (1 + Yq) * [K^ -2s* K%£> u ) + log 

i — Y- i — y_ 

- 2 log (l + Y 1]iVW ) ★ (s * K yU + s) + 2 log — * s * p . v . K vv f x + 2 log — * s 

1 — I I 1 — I_L_ 



u,u 



- log l —^*K l - log (1 - -L) (1 - (4.16) 

where for model I we have 

Vq, = - 2log S* p . v K yU (r^,u) - 2logS(r^ -i/g-u) (4.17) 

+ 2 log Res(ff) * K™' x (u-,u) - 2 log - X + - ^ , 
g x + — x 

and for model II we have 

Vq, = - 2log Si p . v K yU (±rY,u) - 2 log 5(±rf - i/g - u) + (4.18) 

+ 2 log Res(S) * K^ x (u~ ,u)_ - 2 log - 



+ 2 log S-k p . v . K vw * x {-u ,u) + log S vw * x (-u,u) + log S(-u ,u) + 
- log 5i(-tt~,it) - log^i^-tt",^) . 

In the above we use the notations 

/+oo _ , 

dtlog S(u--t)(t-u) K% x (t,v), 
-oo 

s(u) = s(u~) , 

~ x s (u+-) 

the momentum of the magnon is p = i£q{z^) = — Hog — - — fr, and the second argument 

x s\ u —-) 



in all kernels in Q4.16 ) is the Bethe root 



u. 



For a single magnon, the exact Bethe equation in the sl(2) sector can be written in 
the form 

7ri(2n + 1) = i J p — log S s ^rX ( n ; u) +iH, 

where 1Z encodes a correction to the asymptotic Bethe equation represented by the re- 
maining terms. For asymptotic solution we should find 7£ asym = 0, which also serves as a 
non-trivial check of the correctness of the analytic continuation procedure. 

Discarding for large J the exponentially small Yg-functions in the right hand side of 
equation ( |4 . 1 6| ) and subtracting from equation ( [4.16 ) the contribution of the asymptotic 
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Bethe equation, we find 7£ asym for models I and II. We recall that the difference between 
L TBA and J that enters 7£ asym differs between the models; L TBA — J = for model I and 
L T ba — J = 1 for model II. Substituting the asymptotic Y-functions into the respective 
7£ asym , we find that the corresponding quantity perfectly vanishes for both models. 

5. Finite-size corrections to the dispersion: Model I 

For a fixed J and small g, the Yq-functions are small. The finite-size correction term then 
provides the so-called wrapping correction to the energy, or, equivalently, to the scaling 
dimension of a gauge theory operator at weak coupling. The goal of the next two sections 
is to compute the leading wrapping corrections to the asymptotic energy and momentum 
of a single magnon for models I and II. 

When expanding energy and momentum for small g around the asymptotic solutions 
Yq, we obtain the leading order corrections to the energy and momentum 

^o = -^E/d,y Q », P,o = -±±Jdv-^-,Y$(v). (5.1) 

Q=1 J Q=1 J ^ 

We recall that formula ( pM| ) for the energy follows from Liischer's perturbative approach 
and it was shown to successfully describe the four-loop corrections to the anomalous di- 
mension for twist two operators [4, 7]. When the level-matching condition is imposed, the 
functions Yq(v) are symmetric with respect to v — > —v and therefore, Plo vanishes, as it 
should. For a single magnon however, the Yq(v) are no longer symmetric which leads to 
non-trivial values of -Plo- 

5.1 LO wrapping correction to the energy and momentum 

In model I we deal with periodic boundary conditions for all fields (no twist). Using the 
untwisted transfer matrices Tq i, we find that in the limit g — > the lQ-functions associated 
to a single magnon with rapidity u are 

g 2j 4Q 2 (u - Vu 2 + 1) 2 (v - uf 

Y Q( V ) = {y 2 + Q2)J [( Q _ 1)2 + (w _ u) 2][(Q + 1)2 + {y _ n) 2] • ( 5 - 2 ) 

The energy integral ( |5.1|) can be computed by residues for any fixed value of J and the 
first wrapping correction to the asymptotic dispersion relation is then found by taking the 
sum over the bound state number Q. 

In appendix C we present the full result of the leading order wrapping correction 
depending on the rapidity u. Here we present only the result upon substituting the solution 
of the Bethe equations e ipJ = 1. We find for model I the following compact result for J > 4 

El (J) = g 2J £ C m (J)C(2J-2m-l). (5.3) 

m=l 
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Here the coefficients C m (J), written in terms of particie momentum, take the form 

r, its (-l) m+1 r(2J-2m) 2 P ( ■ P\ 2m ~ l , K A , 

CmW = 2 2J - 2m ~ 1 T(J)T(J -2m + 2) tan 4 ( ^ 2 ) X ^ 



x 



(m — 1) sin ^ cos mp — (J — m) cos ^ sin mp 



The above formuia does not hoid for J < 4 because of divergent terms. However, for these 
values of J the wrapping corrections are easily found to be 

£[ (4) = 2C(3)sin 6 |cos 2 |(2 + 4cosp + cos2p) - -£(5) sin 4 |(1 + cosp) (5.5) 
4o(3) = -^C(3)sin 2 ptan 2 | (5.6) 
^o(2) = sm 4 |cosp. (5.7) 

It is worth stressing that for generic p the integral for E^^J) also gives rise to polygamma 
functions. However, these cancel out on solutions of the Bethe equation leaving only a sum 
of (^-functions. 

We conclude this section by mentioning one important subtlety concerning the level- 
matching condition. As we have seen from the single magnon example and as can be verified 
for the case of several particles, without level matching imposed the wrapping correction 
in model I expands at weak coupling as 

E{q = g 2J E i2J) +... = g 2L TBA E (2L TBA ) + _ ^ ( 5 _g) 

since, as we have shown in section 3, in the level-unmatched case L TBA = J. On the other 
hand, with the level matching imposed, one finds that the leading term of order g 2J cancels 
out 15 and the expansion of the wrapping corrections starts from g 2J+A . However, as we have 
seen, for the level-matched case the TBA equations require to set up L TBA = J + 2. Thus, 

El = ^+4^(27+4) +... = g^TBA E {2L TBA ) + _ _ _ ^ ^ 

i.e. in both cases the expansion starts from g 2L TBA^ 

Now we present the leading finite size correction to the momentum. To make the 
relation to Elo more transparent, we introduce 

4o- ( -^V ^- "-^V (5.10 



From the identity 

(u-V^TT) 2 ' v 2 + Q 2 (u - v^TT) 2 ' (u-v^TT) 2 



du Q + " Q + Q =0, (5.11) 



15 The transfer matrix Tq,i is constant as g — > in the level- unmatched case and behaves as g 2 when the 
level matching is imposed. 
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it is then easily shown that 

J&> = -d u S. (5.12) 

It is now straightforward to find the exact expression for the wrapping corrected momentum 
by differentiating the formula ( 8.11 ). Finally, we point out that formula ( |5.12 ) implies the 
following relation between the leading corrections to the energy and momentum 



JPlo = 2 sin ^ sin ~^Lo ~ E lo)- (5-13) 

We would like to note that on solutions of the Bethe equations, the expression for P^o 
generically contains polygamma terms. 

5.2 Quantization of momentum and exact Bethe equation 

To test the equivalence of the two quantization conditions Q3.7| ) , we now need to indepen- 
dently compute the leading finite-size correction to the asymptotic Bethe equations. We 
will do this analytically by applying the method developed in [22, 25]. It is based on finding 
a connection between Yq and the transfer matrix of the Heisenberg spin chain. For a single 
excitation with rapidity u the transfer matrix of the Heisenberg model is 

t M (v) = (M + l)(v-u). (5.14) 

It is readily checked that it satisfies the equation 

t M (v + i)t M (v -i) = t M +i(v)t M -i(v) + t (v + (M + l)i)to(v — (M + l)i). (5.15) 

The corresponding Y-function is then defined as 

n = t M +i(v)tM-i(v) M(M + 2){v -u) 2 

VM(V) t M {v + i)t M {v-i)-t M+l {v)t M -i{v) (M + lf + iv-uf { - > 

which, in turn, satisfies the functional equation 

Vm{v + i)yu{v - i) = [1 + y M +i(v)][l + y M -i(v)]. (5.17) 

Obviously, each um{v) has a double zero at tm = u. Furthermore, vm{v) also enjoys the 
normalization condition 

y M (u-i) = -l. (5.18) 

In terms of the quantities introduced above, the leading order of Yq(v) can be conveniently 
written as 



W = - (t ,2 4 +Q2 )J (V / l+^-") 2 JQ. (5-19) 



2,1 



with 



■ = tQ-i(v) 2 ( z 2 n) 

JQ t (v + (Q + l)i)to(v + (Q-l)i)to(v-(Q-l)i)to(v-(Q + l)i)- K ' ' 
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Moreover, expanding Y^\ vw (v) at weak coupling, we find that the leading term coincides 
with vm{v)- This is all that is needed to repeat the computation of [22, 25] for the case at 
hand. 

Following the same strategy as employed in [21], we employ the exact Bethe equations 
in their hybrid form ( [4.16 ). These equations will allow us to find the leading finite-size 
correction 5u to the asymptotic magnon rapidity u = u(p), the latter given by formulae 
Q3.3| ) and ( |3.4[ ). Expanding ( |4.16| ) around its asymptotic solution, we get 

2vm = Jp + 5K + 0(g 2J+1 ). (5.21) 

Here 57Z is a leading correction of order g 2J to the asymptotic value of 7Z, the latter being 
identically zero. As in the two-particle case [21], the term 57Z is naturally split into three 
parts 

571 = 57l w + 57l {2) + 57l {3) , (5.22) 

where 

^) = -g/d»y g » (Q + 1) " 2 ; (o _ M)2 , (5.23) 

^ <2> - ;E/d»%iw {*><» - -) - g ; ( ;:. u)2 } ■ <*■«> 

The function J-q is defined as 

™ --{{* m - * { ^) - * ( ^) + * ( ^) } ■ <"»> 



The third and most complicated piece is computed in the same way as in [25] by exploiting 
the relation to the Heisenberg model discussed above. Here we only quote the corresponding 
result 



1 °° r°° ( 
*fc(3) = -J2 dvY$ +1 (v) -F Q {v -u) + 

71 Q=l J -° C 1 



(v - u) 2 - Q 2 



(5.26) 



(v - u) ((v - u) 2 + Q 2 ) 
Finally, summing the corrections up, one can show that the resulting expression equals to 



1 00 

5Ti = -y 

2tt 



J2J 



dv 



a 



{v 2 + Q 2 ) 3 dv 



Tq,i(v\u) 
Sq(v,u) 



21 



(5.27) 



which is in agreement with the proposal by [5]. The quantities Tqi(v\u) and So(y,u) are 
given in appendix A. The expression for 57Z can be rewritten as 



571 



1 oo 

— Y 

2vr ^ 



Q=l 



dv 



dYA(v) 2v 

+ J^^YS(v) 



dv 



Q 2 + v 2 



JPlo- 



(5.28) 
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Thus, at leading order the corrected Bethe Ansatz equation takes the form 



2nn 



P + Plo(p)- 



(5.29) 



This provides a strong argument for compatibility of the exact Bethe equations with the 
definition of the exact, i.e. non-asymptotic momentum. Equation (5.29) can be solved 
perturbatively by assuming that the asymptotic value of p is shifted p — )• p + g 2J 5p but 
we will not pursue finding an actual solution any further. Rather, we will continue with 
model II. 



6. Finite-size corrections to the dispersion: Model II 

In this section we will discuss both the leading and the next-to-leading order finite-size 
(wrapping) correction to the energy in model II. 

6.1 LO wrapping correction to the energy and momentum 

In model II the periodicity conditions for bosons are twisted by the total momentum of 
a state; a = p/2. Using the corresponding twisted transfer matrices one finds that for a 
single excitation in the limit g — > the leading term of Yq(v ) is given by 

yo (v)= ^ 16 gdl^ (61) 

Q{ ' (v 2 + Q 2 ) J + 2 (1 + u 2 ) 2 [(Q - l) 2 + {v- u) 2 } [(Q + l) 2 + (v - uf) ' { ' 

The expansion of Yq starts at order g 2J+4 , which also defines the lowest order of the 
corresponding wrapping correction to the Bethe Ansatz energy. Performing a similar com- 
putation of Yq for the multi-particle case, one also finds that in both level-matched and 
level- unmatched cases the expansion of Yq starts from the same order g 2J+4: . This behav- 
ior should be contrasted with the one for Ig-functions in model I based on the untwisted 
transfer matrix. There we found that for the level-unmatched case the corresponding ex- 
pansion starts at lower order g 2J , irrespective of the number of particles being even or odd, 
while for the level-matched case the leading term is of order g 2J+4 . 

In this context it is interesting to recall that in [8] an analytic continuation of the 
leading wrapping correction from even to odd number of particles was proposed. This 
continuation 16 amounts to modifying the (untwisted) transfer matrix Tn % in such a way (by 
multiplying certain eigenvalues of the corresponding monodromy matrix by extra factors 
of i) that in the level-unmatched case for odd number of particles the expansion of the 
wrapping correction starts starts from g 2J+4: , i.e. as for physical states. This modification 
was proposed to describe an unphysical magnon with momentum p = ir. It is easy to see 
that our construction of the twisted transfer matrix agrees and correctly generalizes the 



16 Further subtleties concerning the analytic continuation from even to odd number of particles have been 
discussed in [48, 66], where it was argued that the proposal by [8] corresponds to one (— ) of the two possible 
prescriptions for analytic continuation of the harmonic sums with negative values of indices arising on the 
gauge theory side. 
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proposal [8] for the case of an off-shell theory with any number of excited states and for 
any value of the coupling constant 17 . 

For a generic rapidity it, the resulting wrapping corrections are given in appendix C. 
Here we present the wrapping correction on solutions of the Bethe equations. Similar to 
model I, we find that upon substituting the solution of the Bethe equation u = cot ™, 
the contribution of polygamma- functions in the energy correction completely drops out. 
This again results in a compact form of the wrapping correction, given solely in terms of 
^-functions 



C m { J)((2J - 2m + 1) - 8 



sm 



l p 



m=l 



r(j + |)c(2j + i) 
r(j + 2) ^ 



(6.2) 



Here the coefficients C m {J) written in terms of particle momentum take the form 

Cm(J) - ( lj r(j + 2)r(j-2m + 2) sm ^J* 

msin(m + ^)p — (J + 1) cos | smmp 



(6.3) 



The above formula holds for J > 2 as for J = 2 it contains a divergent term originating 
from C(l)- For this case ££0(2) can be computed directly and the results is 



£&)(2) = 2sin 6 |(2 + cosp)C(3) - \ sin 4 |c(5). 



(6.4) 



It is readily seen that for p = the wrapping correction vanishes. 

Next we present the leading finite size correction to the momentum. From the identity 



d v Y%(v) + a*Y$(v) + 2V q [ J ^ v 2 J y$(v) + 



Au 



4(v - u) 



l + u 2 Q 2 -l + v 2 



Y$(v) = 0, (6.5) 



it is easily shown that 



4 it 



(J + 2)P[ l Q = - d u E* Q - Yf^E^o ~ 4Zi 



where 



(6.6) 



1 /* 



v — u 



+ v 2 — V? 



(6.7) 



The integral X\ can be computed explicitly and it is given in appendix [3,4.2 



17 The twist must necessarily be p-dependent. Indeed, if one admits that in the single magnon case the 
constant twist in [8] is valid not only for p = tt but rather for all values of p, then for the physical state 
with p — one finds a non- vanishing energy correction which is nonsensical because the corresponding state 
being the susy descendant of the BPS vacuum is not renormalized. 
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Using the explicit result for I±, it is now straightforward to find the expression for the 
wrapping corrected momentum by differentiating the wrapping correction to the energy. 
We define the coefficients D m (J) to be 

2 2m - 2J+3 (-l)"T(2J-2m + 2) / . p\2m+4 

Dm{J) = r(j + 3)r(j-2m + 2) l sm 2) 

(J - ra + l)(2m+ 1) sin(n + l)p sin| - 2( J + 2)m sin | cos(n + ^)pj , (6.8) 
and then P^o becomes 



pii _ ( sin l) J ' ~ 



sin 4£ A( J + 2) + sin sin | A{ J + 1) + sin ^£ sin 2 j) 



+2 



4i 2 - J (J + 2)! 2i 1 - J (J + l)! i- 2 - J J! 

£ J D m (J)C(2J-2m + l). (6.9) 



sin I£+i)P s in 3 HA(J-r n - V - 



i-z- J {J- 1)! 



m=0 



Notice that, in contrast to the wrapping correction to the energy, the polygamma functions 
here do not cancel on a solution of the Bethe equation. 

6.2 Quantization of momentum and exact Bethe equation 

Once again, following [22,25], it is possible to establish a connection between the leading 
order of Yq and the transfer matrix of the Heisenberg spin chain. For a single excitation 
with rapidity u we define the following transfer matrix 

t Q (v) = (Q + 1)(Q(Q + 2)+v 2 - u 2 ) . (6.10) 

It is readily checked that it satisfies the equation 

*<j(w + i)H v -i) = *Q+i(«)*Q-i(") + t (v + (Q + l)i)*o(t> - (Q + l)i) • (6.11) 

The associated Y-functions are again defined as 

/ x = *Q+i(")*Q-ifo) /g 12 x 

WJ t (« + i)t (t;-i)-^+i(«)to-i(w)' 

Explicitly, 

[Q2 _ j + ^2 _ U 2][(Q + 2)2 _ ! + v 2 _ u 2] 



I/Q(w) = W + 2) 



[(Q + l) 2 + (w - u) 2 ] [(Q + l) 2 + (v + u) 2 ] • 
These functions satisfy the functional relation 

Vq{v + -*) = [! + I/Q+iOOH 1 + I/Q-i (*>)]■ ( 6 - 13 ) 

Obviously, each uq(v) + 1 has four roots ±uq ± i, that is 

j/ g (±«Q±i) = -l, (6-14) 



where uq = 'u 2 — Q 2 — 2Q. 
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In terms of the quantities introduced above, the leading order of Yq(v) can be conve- 
niently written as 

Y °(v) - 16g2J+4 j Q (a i<\ 

W ~ {v 2 +Q 2 )J+ 2 (1 +U 2)2 > 

with 

• = *Q-i(u) 2 / filfi N 

3Q T{v + {Q + l)i)T{v + {Q-l)i)r{v-{Q-l)i)T{v-{Q + l)i)' K " ' 

where 

t(v) = v — u . (6-17) 

Expanding Y^ v ^-functions at weak coupling, we find Y£^ vw (v) = vm{v)- We would 
like point out that curiously enough yu coincides precisely with the leading term in the 
expansion of Yj^. for the case of two particles with the level-matching condition imposed, 
cf. equation (3.11) in [17]; in particular, each dm has four zeroes at weak coupling. Fur- 
thermore, the transfer matrix tg(«) given by equation ( 6.10| ) which we used to express the 



asymptotic Y-functions for a single and, therefore, unphysical magnon, coincides with the 
two-particle transfer matrix of the Heisenberg model 

tf = (Q + 1)[Q(Q + 2) + (v- Ul )(v - u 2 )] 

with the level-matching condition u 2 = —u\ imposed. 

Using this relationship with the Heisenberg model, one can repeat the computation of 
the leading correction to the asymptotic Bethe equations. This time we find 



1 f 

Q=1 J 



oo 

V — U V — U V + U 

+ 



v° 



.(Q-l) 2 + (v-u) 2 (Q + l) 2 + (v-u) 2 Q 2 -l + v 2 -u 2 . 
The relation with to Plo is 

8K = (J + 2)fgj + 2Xi . (6.18) 

Thus, for model II the leading correction to the asymptotic Bethe equations does not agree 
with the corresponding correction to the momentum, as defined by equation ([O]). This 
result suggests that in a twisted theory the expression for the actual momentum requires a 
modification of expression ( ]3.6| ) in order to agree with the correction to the Bethe ansatz. 
It would be interesting to find its non-perturbative definition. 



6.3 Comparing energies of st(2) and su(2) states 

It turns out that the momentum twisted versions of the st(2) and su(2) transfer matrix are 
related in a very direct manner. In the su(2) sector we will denote the different momentum 
twists by + and — in each factor, writing for a transfer matrix twisted with a = 

±p/2. As discussed before, for s((2) this distinction is not needed as the direction of the 
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twist is inconsequential and we will simply write 1q ^. Expanding the transfer matrices 
at weak coupling, one finds at leading order the following very simple relations 

^ = (q4?) T ^ 2,± <« 9 > 

Now we know that a similar prefactor occurs in the expansion of the Ig-function ( |2.17| ), 

Y Q={gT^) J TlTr ^ + ---, (6.20) 

which means that we can effectively absorb this prefactor in a shift of J depending on 
our choice of twists. In order to indicate this clearly, we write Yq for the momentum 

twisted Y-function in the st(2) sector, and Yq^^ for the Y-function in the su(2) sector, 
where the twists are taken =tp/2 in the left and right factor, respectively. We then find the 
relations 

Y* I(2) (J) = Y* u(2)±T (J) = Y™ {2)++ {J + 1) = Y™ (2) —(J - 1). (6.21) 

Indicating the twists in differently twisted theories by subscript signs, these relations imply 
that the LO wrapping corrections to the energy of the gauge theory operators of the 
schematic form ( , y n DZ J ), (j± T &Z J ), and ( r y±±QZ Jzizl ) all coincide. Here 7 is an element 
of a representation of the orbifold group, appropriate for the field-theoretic realization of 
the twist, and the integer n labels the n-th twisted sector in the orbifold theory. Finally, 
the twists ++ and also affect the su(2) Bethe equations, cf. equation ( 2.22j ). Now we 



immediately see that these twists can be effectively compensated by shifting J — > J ± 1, 
leading to the one and the same Bethe equation e ipJ = 1. Therefore, the theories obtained 
by doing the various twists are directly related by a simple shift in J in the su(2) sector. 



Finally we would like to note that there exists a curiously simple relationship not only 
veen 
Namely, 



between the leading order but also the full twisted transfer matrices ^ and T^y . 



T s Q T +T Sf y = T Sf- ( 6 - 22 ) 

An analogous relation holds in the Konishi case, where it was observed that for level 
matched states T^y and coincide [17]. Of course this agrees with the above relation, 

since for level matched states the twist is trivial 18 . 

6.4 NLO wrapping correction to the energy 

To describe the next-to-leading order wrapping correction all quantities in fl3.1|) must be 



carefully expanded to two orders higher in g. In what follows we will discuss the expansion 
of the different terms separately and give their contribution to the energy. 

First we consider the term in the integral that is not related to the Y-function 

i + V , V l~^L +0{g A ). (6.23) 



dv a {v 2 + Q 2 f 



18 For a single magnon the level matched equivalence is completely trivial. 
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This obviously gives rise to the contribution 

,2 ra 



<o - V E / dv YqW- (6-24) 



This term can be written as the sum of two terms 

E§lo(J) = (J^yi " 2i&(J + 1) (6.25) 

The new integral I2 can be expressed in terms of the previously computed wrapping cor- 
rections. The explicit expressions in this integral are presented in section S.4.2; . 

Next we focus on the term \og(l+YQ(v\u)). This term depends both implicitly (through 
the Bethe root u = u + g 2 5u) and explicitly on g. Expanding this gives 

log(l + Y Q (v\u)) = Y^°(v\u) + Y^\v\u) + Y^°{v\u) + 0(g 2J+s ), (6.26) 

where Yq°(v\u) is of order g 2J+4 and Yq° ,2 (v\u), Yq LO (v\u) are of order g 2J+ ®. Here 
Yq (v\u) is obtained by expanding the Bethe root u = u + g 2 5u and is expressed as 

Y$ LO (v\u) = g 2 d u Y^°{v\u)5u. (6.27) 

The resulting contribution to the NLO wrapping correction from this piece is consequently 
given by 



E^l {J)=g 2 smp^-E^ (J), (6.28) 



where we used (3.3) to express 5u = smp. 

Subsequently, we turn our attention to Yq (v\u). It is useful to split this function 
into two parts 



4v 2 2vu 4 J (v 2 - Q 2 



(Q 2 + v 2 ) 2 {u 2 + 1) (Q 2 + v 2 ) (u 2 + 1) (Q 2 + V 2 ) (Q2 + w 2)2 
4 4(Q 2 -1) 1 4 



v 



(Q 2 + v 2 ) (1 - Q 2 - v 2 + u 2 ) Q 2 + v 2 (u 2 + lf 



YX(v), (6.29) 



=g 2 y 2 ; 2 + 1 2 2 ~j Y&v). (6.30) 

The contribution E^ Q of the first term Y^° can be related, by some partial integrations 
to integrals that were encountered before. In this respect it is computed in a straightforward 
manner. The relevant integrals are collected in appendix D. The second term comes from 
expanding the dressing phase, as presented in appendix A. Note that we here used harmonic 
numbers H n rather than the digamma function ip(n). They are related via H n = 7 + 
i[)\°)(n + l). Appendix D is devoted to computing the corresponding contribution which 
we denote as E^q. 
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Finally, we consider the term involving the dispersion relation £{p) from equation (|3.1|), 
Since the momentum also receives a wrapping correction, i.e. 

p^p + g 2J+A 5p, (6.31) 

the asymptotic energy £{p) also gets corrected 



8( p ) = yi + 4 5 2 sin 2^±^ = / : + 4^2 gin 2 | + g 2 J+6 gin p Jp + ^2 J+8 j _ ( 6-32 ) 

The TBA approach relates <5p to the correction to the asymptotic Bethe equation 51Z as 

5p = --^-. (6.33) 

L TBA 

Now it is time to recall that for model II the correction 5p is not proportional to P^ . On 
top of this, from our TBA considerations in section 4 we found that Ltba = J + 1, which 



results in the following expression for 5p, cf. formula ( 6.1 



Before continuing with our discussion of the next-to-leading order correction to the energy, 
let us mention that there is a consistency check on our result. Each of the above discussed 
intermediate terms contain polygamma functions. However, just as for the first wrapping 
correction, we can expect them to cancel in the final result. It turns out that the contri- 
bution for polygamma functions indeed exactly cancel out for our rather non-trivial result 
for 5p. 

Collected results 

We have explicitly computed the next-to-leading order wrapping corrections for states 
with J = 2,3,4 and for one state with J = 6. We found that in each case the result is 
a combination of £-hmctions and multiple (double) zeta- functions (see appendix D). We 
will present the results with all double zeta values expressed in terms of usual C-functions. 
Below we summarize our results. 



For J = 2, the Bethe equations admit only the solutions p = 0, it, i.e. u = oo and u = 0. 
The contribution for p = vanishes, and the NLO wrapping correction for p = it is found 
to be 

" -4C(3) - jk(3) 2 - 5 C (5) + ^C(7). (6.35) 



glO ^ 2 ' 

J = 3 

For J = 3 the Bethe roots are u = db-4=. Both rapidities yield the same wrapping correction 

V3 

= 512 C(3) " 128 C(5) " 512 C(3)C(5) ~ + ^ C(9) - (6 - 36) 
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J = 4 

For J = 4 there are three Bethe roots u = 0, u = ±1. For w = 0we find 

^# =C(3) 2 - 4C(5) 2 + C(3) f ^(5) " tC(7)1 (6-37) 



^ ' 32 



while for u = ±1 the result is 



- £<(6)-7C(7)-^C(9) + ^C(ll), 



9 >/r , 119 2037 H55 A/11N 



J = 6 

Finally we computed the NLO wrapping correction for J = 6, u = and found 

_ 5 C(7) _ 7 c(7)2 _ c(5) / 3 c(?) + 15 c(9) \ _ W c(g) _ 957 c(11)+ 



gl6 2 M y 512 sw sw ^32 64 sv 7 8 sw 128 

(6.39) 

+ c(3) (|c( 5 ) + > + 2c( 9 ,-Sc(n,)-^C( 1 3, + ^c( 1 ,. 

One can see from these explicit results that the term of maximum transcendentality is 
always of degree 2 J + 3. In case of the leading order wrapping correction this was 2 J + 1. 

In the above, we have chosen to express the double zeta-functions in terms of usual ('s 
and products thereof, but we can also attempt to go the other way around and express the 
NLO correction purely in terms of double zeta-functions; this turns out to be possible. As 
a matter of fact, all the wrapping corrections are of the following form 

E% LO = g 2J+6 ^ A v -C(i,j). (6-40) 

i>j 

For instance, for J = 2 this formula explicitly reads 

£nlo(2) = - 4C(2, 1) - 10(3C(4, 1) + C(3, 2)) - |(16C(5, 1) + 9C(4, 2))+ 
1 05 

+ — (20C(6, 1) + 10C(5, 2) + 3C(4, 3)). (6.41) 
Since the leading order wrapping correction was of the form 

£lo=5 2J+4 E^«, ( 6 - 42 ) 

i 

the general structure that is expected to be found for higher order wrapping corrections 
will be ordered sums of multiple zeta values. 
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We close this section with the following curious observation. As we have seen, the 
relation between 5p and Plo is given by equation ( 6.34| ). It is interesting to ask what 



we would find for the NLO energy correction if we would assume that bp and Plo would 
be related in the standard fashion: 5p = —Plo- It turns out that the contribution of 
ifj- functions also cancels out completely in the expression for Pnlo! The final result differs, 
however, from PjjLO' om y i n terms linear in ("-functions. For example, for J = 4, u = 1, 
the difference is 2 {340 — hr fact, we have found that the following linear combination 

Pl 1 o+21 1 (J) (6.43) 

is "^-free" and, therefore, being added with arbitrary coefficient a to our original expression 
for 5p does not spoil the "^-free" property: 

dp = -P?o + 7^( P lo + 2Xi (J)) • (6.44) 

When u = both choices for 5p give rise to the same NLO correction to the energy. The 
implications of these observations are not clear to us. 

7. Finite-size corrections in /3-deformed and physical orbifold theories 

In this section we apply the generalized Liischer formulae to compute the leading finite-size 
corrections to the energies of one- and two-particle excited states from the su(2) and sl(2) 
sectors in /3-deformed theories, and to two-particle states in the sl(2) sector of a generic 



Zg-orbifold. The Yq (u)functions needed to compute Elo according to equation (l5.l|) are 



constructed by using the twisted transfer matrices and the general formula ( p. 17 ). The 



rest of our considerations follows the same logic as for models I and II. We also make a 
comparison between our results for /3-deformed and orbifold theories for a single magnon. 

7.1 One-particle /3-deformed states 

For single particle states, we find the following energy corrections for the su(2) and sl(2) 
sectors. 

su(2) sector 

According to our discussion in section 2, a single magnon in the su(2) sector of /3-deformed 
theory has the following Bethe equation 

1 = e i<p-W)J ^ p = 27T f3 + ^, n = 0,...[|l. (7.1) 

For the level-matched case n = so that a physical magnon has p = 2tt(3. Its asymptotic 
energy is 



E = J + y 1 + Ag 2 sin 2 vr/3. (7.2) 
On the gauge theory side an operator corresponding to the physical magnon is tr($Z J ). 
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Let us denote the transfer matrix obtained by twisting the bosonic eigenvalues (in the 
auxiliary space) with a as T(v\u; a). Implementing the twist describing /3-deformed theory 
results in the following Y-functions 



■2.1 



Yq(v) 



T s ^ i (v\u ] p)T s ^ r (v\u;(2J-l)P) 



(Q 2 + v 



2\J 



S (u,v) 



(7.3) 



where the coefficient So(u, v) is given in appendix A. The leading term in the weak-coupling 
expansion for the twisted transfer matrix in the su(2) sector is explicitly given by 



(v\u;a) 



VQ 2 



Qe 



iiza. / u-\-i 



^(Qz + ( u -v + iy) 



(Q-l){(Q + lY + ( v -uf) 

(i(Q - 1) + v - u)(i - u) (i(l - Q) - v + u){i - u) 

Qe-™ a (i(Q + 1) - v + u) (Q + l)(i(l -Q) + v-u) 



-+ 



^ / V u— % 



(7.4) 



This transfer matrix starts at order g which results in the fact that our Y-function, and 
corresponding corrections to the energy and momentum, will be of order g 2J+2 . 

It is now straightforward to compute the first wrapping correction to the energy. We 
find it to be a sum of (^-functions 



El/8 ( T\ — r, 2J + 2 



where the coefficients are given by the following expressions 



B^J) 



■ sin(Jvr/3) sin^vr/?) sin(( J - 2)tt/3), 



(7.6) 



B^J) = sin(Jvr/3) sin^ n (7r/3) [n - 1) sin((J - 2n)vr/3) - Jcos(tt/3) sin((J - 2n + l)vr/3) 

The leading finite-size correction to the energy found above perfectly agrees with the field- 
theoretic result obtained in [44,45]. As discussed in the introduction, this is also in agree- 
ment with the results of [46]. 

Special attention has to be paid to the particular case of rational /3, (3 = j with 
n = 0, 1, . . . ,J — 1. For these values of the deformation parameter, the naive lowest order 
in the weak-coupling expansion of the transfer matrix Tq^\v\u; (2 J — l)/3) vanishes and 
the expansion really starts two powers of g higher. Bearing in mind that the rapidity u 
is g-dependent via (|3.3| ) with momentum p = 2irf3, it is readily checked that the functions 
Yq for these special values of f3 are absolutely the same as for generic f3, up to a shift 
J — >• J + 1, so that we find 



E\ 



fi=n/J 
LO 



(J) = EUJ + 1 



(7.7) 



In particular this implies that for these special values of f3 the expansion of the Yg-functions 
starts at order g 2J+A . 
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sl(2) sector 

In the sl(2) sector the single magnon Bethe equations are e LpJ = 1 so that 

27T77- 

P= — , n = 0,l,... J -I. (7.8) 

The physical state must be level-matched leaving only the solution p = 0. In the sl(2) sector 
of /3-deformed theory the left and right transfer matrices Tq^ are twisted with a = and 
a = 2J/3, respectively. Computing then the corresponding twisted Y^-functions, we find 
that they all vanish for p = 0, resulting in a vanishing wrapping correcting to the energy. 

7.2 Two-particle /3-deformed states 

The twist should correctly account for /3-deformation of N = 4 SYM for any number of 
excitations; in this section we will briefly investigate some two-particle states. In particular, 
we reproduce the wrapping energy correction found in [45, 47] for the Konishi state in the 
su(2) sector of /3-deformed theory, and investigate two-particle states in the st(2) sector for 
which no field-theoretic results are known so far. 

su(2) sector 

The perturbative computation of the anomalous dimension of /3-deformed su(2) Konishi- 
like states [45] agrees perfectly with results coming from the asymptotic Bethe equations 
up to three loops. At four loops, there is a discrepancy arising due to wrapping effects, 
which should accounted for via energy corrections computed from the Liischer formulae. 
In [47] a modified S-matrix was used to obtain an expression for the wrapping correction 
to the energy, which is in agreement with the explicit perturbative results. 

The modification of the S-matrix proposed in [47] amounts to giving the matrix ele- 
ments of the S-matrix for fermionic states a /3-dependent phase: 



This is obviously different from the twisting procedure we have proposed above, and there- 
fore we cannot a priori expect these two different approaches to agree. However, the su(2) 
Konishi state is special, because it has particularly symmetric quantum numbers; J = M. 
This means that for such a state our twist amounts to a left twist by M and a right twist 
by 2 J — M = M. For such a symmetric twist it is not hard to see that it can effectively 
be absorbed in the unmodified S-matrices, exactly in the fashion indicated above. This 
means, in particular, that for the Konishi state the leading wrapping correction computed 
from the twisted Y- functions must naturally agree with the result of [47]: 

81(1-3A) 2 (1 + A) 4 l 



^Konishi — B 



-54(1 + A) d (-5 + 3A)C(3) - 360(1 + A) 2 C(5) + 



(1 + 3A) 2 



where A = yJEEI|<1iIze2 . Of course, we have also verified this agreement by explicit calcu- 
lation based on the twisted transfer matrix, both analytically for certain nice values of /3, 
and numerically for various random values of /3. 
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We would like to emphasize that only in the special case of J = M it is possible 
to effectively absorb the twist in the modification of the S-matrix in a natural way, and, 
therefore, these different approaches will only agree in this special case. 

st(2) sector 

Now we apply our approach to analyze the wrapping corrections in the s[(2) sector of /3- 
deformed theory. We restrict ourselves to two-particle states with angular momentum J, 
corresponding to the gauge-theory operators schematically of the form tv(D 2 Z J ). 

To obtain the corresponding lQ-functions, cf. our earlier discussion in section 2, we 
need to twist the right transfer matrix with a factor of 2Jirf3 and to keep the left trans- 
fer matrix untwisted, regardless of the number of excitations. However, since the Bethe 
equations in the s[(2) sector are not modified under /3-deformation ( [2.19|) , the standard 
level matching condition remains p = 0. With this twist and the level matching condi- 
tion ui = — U2 = u taken into account, the left transfer matrix for two particles take the 
following form in the weak coupling limit: 

s i(2)e = 8g 2 Q (Q 2 - 1 + v 2 - u 2 ) 

G' 1 (1 + u 2 ) (Q 2 + v 2 ) {{Q - l) 2 - 2i(Q - l)v -v 2 + u 2 ) ' 

which for J = 2 is of course nothing else but the transfer matrix for the usual Konishi 
state. On the other hand, due to the twist, the right transfer matrix starts with a non- 
trivial zeroth order term 



^1(2) r _ 4Qsin 2 (vrJ/3) (l - Q 2 + 3v 2 - 3u 2 ) 

" * ( u + i (Q_l))2_ u 2 > 



This immediately implies that the wrapping correction will show up at an order g 2 lower, 
except for the special values 19 (3 = m/J, m = 0, 1, . . . , J. Putting these results together in 
equation ( 2.17 ), yields the lQ-function 



YQ = sm mnnr- ' 

where 

f$ = (Q±l) 2 + (vTu) 2 . 

Thus, we are now able to compute the wrapping corrections to the energy for two-particle 
s[(2) states with an arbitrary twist. Evaluating the energy integral for arbitrary u, gives 
a string of polygamma functions in addition to the expected ("-functions. However, substi- 
tuting the actual solution of the Bethe equations, these polygamma contributions cancel 
out in a highly non-trivial manner. 

Here we explicitly present the wrapping corrections to the energy for a few values of 
J where the rapidities take particularly simple values; since in the small coupling limit 
uj, n = cot jxj, where n = 1, . . . , |_^f-J > we restrict ourselves to J = 2, 3 and 5. The results 
for the corresponding wrapping corrections are summarized in Table 0. 



19 For these values the twist effectively disappears, i.e. we return to the undeformed theory. 



- 39 - 



J 



n 



u 



1/V3 
1 


1/a/3 




2J+2 



sin 2 (Jvr/3)) 



-n-K( 3 ) + |C(5) 

-6 + f C(3) - f C(5) 
~ik " T^C(3) + ^C(5) - iiC(7) + |C(9) 
" IC(3) + 1C(5) + £C(7) 
-f ~ jCg) + fC(5) ~ jC(7) - §C(9) 



Table 3: The leading finite-size correction to the energy for certain two-particle states from the 
sl(2) sector in /3-deformed theory with /3 ^ j, m e Z. For J = 2 the leading correction arises at 
three loops and for J = 3 at four. Note that compared to the other two J = 5 states, the term of 
maximum transcendentality for J = 5,u = 1/V3 is missing. 



7.3 Physical two-particle sl(2) orbifold states 

While the Bethe equations for a two-particle state in the sl(2) sector are not modified by 
the particular orbifolding corresponding to the twist ( |2.5| ), the wrapping corrections to the 
energy are. To consider two-particle states on a Zg-orbifold, we should give both sl(2) 
transfer matrices the same twist; one of the S-th roots of unity. The explicit expression 
for the lowest order term of the twisted transfer matrix has already been given in equation 



(7.10), where now all we need to do is replace 2J-7T/3 by 2irm/S, where m = 0, . . . , S — 1, 



and m = of course gives the usual sl(2) descendant of the Konishi state. Using formula 



( 7.10 ), the leading finite-size correction to the energy of any desired state can be readily 



computed, here we present the explicit results for some simple states in Table |I[ 
7.4 Comparing /3-deformed and off-shell orbifold magnons 



According to equation (7.1), a single su(2) magnon in /3-deformed theories has momentum 
p = 2-7T/3 + 22p. Clearly the case of /3 rational, j3 = j, ra 6 Z, is special, because now 
p = 27r ("^+") ; j_ e momentum is essentially the same as of an off-shell st(2) magnon in 
orbifold theory. This suggests the existence of a certain relation between the su(2) magnons 
of /3-deformed theories with f3 = *j and sl(2) magnons of orbifold theory; their energies 
must be equal. We could also come to a similar conclusion by inspecting table |l|. Indeed, 
for f3 = J the right twist in the su(2) sector is (2J - M)vr/3 = 2vrm - vr/3M = -vr/3M. It 
is clear now that if for model II the total momentum is chosen to be p = 2-7T/3M, then the 
twist in the su(2) sector of /3-deformed theory just coincides with the one in model II. 

In fact, the equality between the energy of an su(2) magnon in /3-deformed theory with 
a rational (3 and the energy of the sl(2) magnon in the off-shell theory (with /3 = 0) has 
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J 


S 


u 


E orh (J)/(g 2J sm^^m/S) 


2 
3 
3 
5 
5 
5 


1 
1 
2 
1 
2 
3 


1/V3 
1 


V3 
1/a/3 



l 

3 

-n-K(3) 

-l + f log(2)-iC(3) 
-BT6 " wC(3) + ikC(5) - §C(7) 

-i-|C(3) 
-f + f C(3) - |C(5) - 



Table 4: The leading finite-size correction to the energy for certain two-particle states from the 
sl(2) sector for a Zs-orbifold. The leading order correction enters at order g 2J . It is interesting 
to note that compared to the other two J = 5 states, the term of maximum transccndcntality for 
J = 5, it = l/\/3 is missing. 



been already pointed out in [48]. There it was argued to arise due to the symmetry of 
the spectrum under the shift /3 — > j3 + j. In particular, a physical magnon in /3-deformed 
theory with (3 = 1/2 was shown to have the same leading correction to the energy as a 
physical magnon with p = tt in the undeformed theory [48,66,67]. The formulae for the 
energy corrections we have obtained allow us to give a proof of the corresponding statement 
for arbitrary J. Indeed, comparing formula ( |7.5| ) taking for ft = j against the wrapping 
correction (|6.2j ) in model II for p = one finds 

E (3 (J + l) = E ll {J). (7.11) 

Actually, the observation made for model II extends 20 to the level of Y-functions. 
Namely, we have that 

cos(tt^) - cos(tt(2 J + l)/3) v-u „ 

sin 3 (7r/3) Q2-l + vi- U 2 Q [J> - 

We see that for (3 = j only the first term on the right hand side of the last expression 
survives. This also automatically implies that P^o f° r @ = 7 exactly coincides with P^o- 

20 One can try to compare /3-deformed theory against model I. Here one does not expect any obvious 
relation, but we still find a certain structural similarity of the leading order corrections to the energy up to 
normalization and the piece of maximum transcendentality 

E?(J - 1) = -64cos 4 E\J) - 8,-sin 4 (^) C (2J - 3). (7.12) 
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Notice that the second term in equation ( 7.13 ) coincides with the difference between 51Z 
and P£ Q . 

Finally, we can speculate on the relevance of the NLO result for model II in the 
context of /3-deformed theories. If for model II the equivalence between the sl(2) and the 
5u(2) sectors continues to hold beyond the leading order, it is probable that the results 
for -EJjlo we f° un d in section 6 also represent the NLO energy corrections for su(2) states 
in /3-deformed theory with (3 = j. Coincidence of P^o w hh P^o can be considered as 
serious evidence in favor of this. However, since currently we do not know what 51Z is in 
/3-deformed theory and how it is related to P^q, we can not directly verify the equality of 
the NLO energy corrections. 

We could also speculate that in /3-deformed theory the situation might be slightly 
better, in the sense that possibly 5p = P^q, in line with our curious observation at the 
end of section 6. In that case the corresponding energy correction could be easily obtained 
from our results for E^ LO and it would differ from the later by terms linear in ^-functions. 
However, for the moment this is just a bold speculation. 



Acknowledgements 

We are grateful to Sergey Frolov for numerous important discussions, and to Zoltan Bajnok, 
Christoph Sieg and Matteo Beccaria for valuable comments on the manuscript. G.A. 
acknowledges support by the Netherlands Organization for Scientific Research (NWO) 
under the VICI grant 680-47-602. The work by M.L. and S.T. is a part of the ERC 
Advanced Grant research programme No. 246974, "Super symmetry: a window to non- 
perturbative physics". 



8. Appendices 



8.1 Appendix A: Twisted transfer matrix 

The eigenvalue of the twisted transfer matrix for an anti-symmetric bound state represen- 
tation with the bound state number Q is given by the following formula, generalizing the 
result of [61] 

K 11 



yi-x^ y x 



i=l 



. JT Vi~ x 



yi-x-<- V x~ 



i=l 

Q-1K 11 



x ++ i _ y ._i_«a 



n[ 

i=i 



(x -x i )(l-x xj) x + 



(x+-xr)(l-x+xt) *" 



X /x" 1 



fc=l 1=1 

/ / 1 1 %-x+ V x_ 



x+^-^-J,- 



if 1 



i=l 



k=0 i=l 

X-iii 

X < 1J 

i=l 



x+ " Vi 9 



J~\\+(v,Ui,k)+ Y[ \-{v,Ui, k) 
1=1 

n 



x ' —x- / x„- 



I T -I; \/ x; 



2ifc 
9 



D-U; + i(Q-l) 



c+ _ 1 i(2fc-l) 



UJj — X T — 



1 i i(2fc+l) 



■+e 



II . 



J__^+ 1 | 2ifc 



Vii X T — ■ 

= 1 Wi 



-^" In _i_ l . i(2fc+3) 



, 1 i(2fc+l) 



> . 



Here the twist e ia enters only the last line. Eigenvalues are parametrized by solutions of 
the auxiliary Bethe equations: 



n 



m - Xj 

Vk - x\ 




1=1 

w k ~ Vi ~ k + I 
t = [ w k - Vi - £ - | 



Wi ~ yk ~ m ~ 9 



— + - 

Vk 9 



t=l,t^k K 1 9 



(8.2) 



2in ir "V- - j- 

2i ' 



In the formulae above the variable 

v = x + + 



1 



% 1 2 

_Q = x - + — + -Q 

9 x g 

,± - 



takes values in the mirror theory rapidity plane, i.e. = x(v ± ^Q) where x(v) is the 
mirror theory x-function. As was mentioned above, Uj take values in string theory u-plane, 
and therefore x^ = x s (uj ± |) where x s (u) is the string theory x- function. These two 
functions are given by 

1 , 



x{u) 



[u — l 



yjA-u 2 ), 



u , / 4 , 

X.(u) = -(l + 4/l-- 5 ). 



(8.3) 



Finally, the quantities A± are 



X±(v,Ui,k) = 



1 



(x i x + — l)(x + — .t^) 2ifc a;+ ( a; i + x t) 



{Xi ~ x+)(x+xf - 1) - z+)(a;+a;t - 1) 



jg+( g< -a;t) / 
(xr-x+)(x+x+-l)\l 



i(2k-Q) 



(8.4) 



-43- 



For a single string theory particle with rapidity u the twisted sl(2) transfer matrix is 
of the form 



T* Q T(v\u)=l 



(■« 



T)(l 



; x ) (l - 2:^2;+) ir 



2 cos a 



Q(w - u) 



(u-«)-i(Q-l) 



+ A, 
(8.5) 



where 



A=(Q-1) 1 



The matrix Tf'fVl 

v) is normalized as 



(x + — xf)(X-y X + — 1) 

(x+ — 2:^) (2^2;+ — 1) 



-0(0 - 1) x+ ^ Xl 



which is immediately clear form its expression above. Here the star indicates analytic 
continuation to the kinematic region of string theory. 

The S-matrix in the string- mirror region 5^,^ is found in [68] (see also [5]) and it has 
the following weak-coupling expansion 



where 



and 



S sl(2)( u ' v ) = S Q (u,v)+g S 2 (u,v) + ... , 



r< 1 , [{v-uf + {Q + lf} [Q-l + Hv-u))} 
So ^ V >- (u-i)*[Q-l-i(v-u)] 



(8.6) 



S 2 (u,v) = -S (v,u) 
S (v,u) 



2[2Q(u + (it + i)(v 2 + Q 2 + 2v(u - i))] 
(« 2 + Q 2 )(l + u 2 ) 



+ 



(8.7) 



1 + u 2 



Q + iv\ ( Q + iv \ / lj — iv . 

4 7 + v 1 + -Z-z — + i> 1 - -t-z — ) + ip [ 1 + — ) + ip ( 1 - 



Q — iv 



These expressions are enough to build up the two leading terms in the weak-coupling 
expansion of the asymptotic function Yq. 

8.2 Appendix B: Asymptotic critical points and TBA 

As mentioned in the main text, a single magnon shows a critical behaviour quite similar 
to what has been observed in the two-particle case in [17]. However, we would like to 
stress that the discussion of criticality is restricted to the asymptotic solution; the data 
recently obtained for the two-particle case [24] is not sufficient to conclusively establish the 
existence of critical points for an exact solution. 

For model I we find critical behaviour that arises from the movement of zeroes of Y + — 1 
and poles of Y M i w and Y M i vw in the complex plane. These zeroes and poles run into the 
integration contour at certain values of g and further cross it, yielding extra contributions 
in the TBA equations above these what are called critical values of g. 
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We indicated before that zeroes of 1 — Y + and poles of Y M \ W and Y M \ VW are related 
by shifts of i/g; owing to this we can treat their critical behaviour in a unified manner. 
For concreteness however, let us first consider zeroes of Y + — 1. The function Y + has the 
interval (—2, 2) as its associated canonical integration contour, and as is plotted in Figure || 
the zeroes of Y + — 1 cross this contour, and do so at a particular value of g denoted g® r . As 




Re(r) 



Figure 6: The behaviour of zeroes of Y + — 1 as g increases, for p taking the values ^qqqq , ^f- , 



7T 7T 

3 ' 5 



Yq, from left to right. For negative values of momenta the picture should naturally be 
mirrored with respect to the imaginary axis. 



these zeroes do so from either side of the real line they drag the integration contour along 
with them. In order to maintain the canonical integration contour in the TBA equations, 
one needs to pick up a contribution from each of these zeroes as they move back into the 
complex plane, as illustrated in Figure [?]. For the poles of Y M \ W and Y M \ VW the story is 



® 



Figure 7: The zeroes of Y+ crossing the real line giving extra contributions as the integration 
contour is returned to the real line. 

nearly identical as they show exactly the same behaviour, except that the crossing happens 
at higher values in g. The critical values as a function of the particle momentum p have 
been plotted for the first few M in Figure ||. The contributions these zeroes give in the 
TBA above the critical values can straightforwardly be read of from the ground state TBA 
equations taking into account the orientation of the integration contour. 
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Figure 8: The critical values g„ as p varies from to ir. 



In model II the story is different as 1—1+ has no zeroes, but now the zeroes vm 
show criticality. As discussed in the main text, we have found exact expressions for the 
asymptotic zeroes, allowing us to understand and treat the asymptotic criticality exactly. 
Concisely the story is as follows; a zero rjy/ can be imaginary but between the real line 
and — i/g at small coupling for a given state (momentum). When this is the case, the zero 
should be accounted in the excited TBA equations, but as g is increased the zero moves 
away from the real line and crosses either the real line or —i/g. The value of g for which 
this happens can also be called a critical value, where now beyond the critical value a zero 
tm no longer contributes to the TBA equations. This story can be readily verified by 
making use of the explicit expression fl4.5|) . Please note that in this twisted model, for a 
carefully chosen state these critical values can be arbitrarily small. 



8.3 Appendix C: Wrapping corrections 

In this section we present the LO corrections to the energy and momentum in full generality 
for both model I and II. 



8.3.1 Model I 

We write t/;( m ) for the m-th polygamma function, ^( m ) = - ^J^ , where tp 
V (z)/T(z). Let us first define 

(u - V^Tl) 2 



2tt 



Then, we introduce the following quantities 
A(J) = J 



B m (J) = 



T(J-m)T (j-m+i) 
r(J- 2(m- 1)) 



to — 1 + i( J — m)u in — 1 — i{ J — m)u 



(1 + iuf 



{l-iuf 



(8.9) 
(8.10) 
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The full leading order wrapping correction to the energy for J > 4 is given by 



5 2J 2- 7 (J-l)!(u 2 + l) J ~ 1 I 4J w 2 + l u 2 + l ; J 

-Tj^V E S m (J)C(2J-2m-l). (8.11) 



Formula (8.11) was checked analytically up to J = 8 and numerically up to J = 200. It is 
also readily seen that the energy correction is real. For J = 3,4 formula ( |S.ll ) is not well 
defined due to the divergence £(1) = oo. This can be resolved by taking the appropriate 
values of J, n and substituting £(1) — > in ( fe-liP - After this substitution, formula (|8.11| ) 
holds. Explicitly, for J = 3 one gets 



S (3) _tt [(iu + l) 3 - (iu - l) 3 ] f A(3) iu 
~f~ ~ 16 (u 2 + if \~12~ + u 2 + l 

7t (1-3u 2 ) 7£ + 3m 2 (1 + u 2 )C(3) 

2 (u 2 + l) 3 



A(2) 



u 2 + 1 



.12) 



For J = 2 one finds (A(0) = 0) 



<f(2) _ 7r[(m + l) 2 -(m-l) 2 ] f ,4(2) 
~~ 4(u 2 + 1) l - ^ - 



u 2 + 1 



7r (1 - 47 £ )m 2 - 1 
2 (u 2 + l) 2 



.13) 



Finally, for J = 1 the sum over M does not converge and we find that the wrapping 
correction diverges. Nevertheless, expression ( |S.ll ) can be continued to J = 1 and it gives 
a finite result 



7T 

12 



4^ 2 
(u 2 + if 



3^ 



+ 3ip 



(8.14) 



So far the formulae for the energy correction were obtained for arbitrary values of 
particle rapidity u, i.e. not necessarily satisfying the quantization condition. If we impose 
the quantization condition (the Bethe equations), then at leading order in g the rescaled 
rapidity is u = cot with n integer. For these values of u the following relation holds 



(iu + 1) J - (iu - l) J = 55-7(1 - e 2 ™) = . (8.15) 



Therefore, the contribution of the polygamma- functions to the energy correction completely 
drops out on solutions of the Bethe equation, resulting in formula (5.3). 



8.3.2 Model II 
Energy 

To present the value of the energy integral for model II, we define the following 

quantities 



A(J) = 



B m (J) = 



r(J-m + l)r(j-m+|) 
r(J- 2(m- 1)) 



m + i(J — to + l)u m — i(J — to + \)u 



(l + iu) 2 



[1-iuf 



(8.16) 
(8.17) 



Then we find for J > 2 



,2J+4 



2 J J!(u 2 + l) 



2(J+l) u 2 + 



u 2 + l 



(8.18) 



(1 + u 2 ) 



WJ+ 1 

2 ^F(7TI)! E ^(WJ-am + i) 

v v 7 m— 1 



8 r(j + |)c(2j + i) 
(l + u 2 ) 2 r(j + 2) 0F ' 



Formula (|8 . 1 8|) was checked analytically up to J = 8 and numerically up to J = 300. 



For J = 2 formula ( |B.ll ) is not well defined, again because of the divergence £(1) = oo. 
However, for J = 2, the energy correction can be easily evaluated directly 



^ I (2)^(m-l) 2 -(iu + l) 2 J J 4(3) 2iuA(2) 4A(1) 1 (6u 2 + 2) C(3) 



(u 2 + iy 



u 2 + 1 u -\- 1 



5C(5) 



(u 2 + iy 



2{u 2 + iy 



Finally, for J = 1 the sum over M does not converge and, therefore, the leading wrapping 
correction diverges. 



Momentum 

Define 



D m (J) =i 



then 



2 2m " 2J+2 r(2 J - 2to + 2) 

r(j + 3)r(j-2TO + 2) 

2(J+2)to ( 1 



{J -771+ 1)(2TO+ 1) 

1 + u 2 

1 



\2m+l 



P, 



LO 



(1 + u 2 ) 2 \(l + iu) 



£ D m (JX(2J-2m+l)- 



(1 - iu) 



2m+l 



{l+iu) 2m + 2 (l-iu) 2m + 2 



(8.19) 



+ i 



(iu - 1) J - (iu + 1) J 
2- / + 2 (J + 2)!(u 2 + l) J+1 
.(m-l) J+2 -(iu + l) J + 2 



A(J + 2) + % 



(iu + l)(iu - 1) J - (iu - l)(iu + 1)' 



2 J + 1 (J+l)!(u 2 + l) 



J+2 



A(J + 1) 



- ,(»» + i)(»»-i) J+2 -(«»-i)(»; + i) J+2 A(J _ 1} . 



2 J J! (u 2 + 1)' 7+3 ' 2^- 1 (J-l)!(w 2 + l) J+4 

On the solution of the Bethe equations the polygamma functions do not drop out 
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8.4 Appendix D: Details on the NLO computation 

In this appendix we have gathered some details on the explicit computation of the NLO 
correction to the energy. 

8.4.1 The Harmonic number integral 

Special attention has to be paid to the contribution Y^ LO due to the appearance of har- 
monic numbers. Let us first set up the tools to handle sums and integrals of these numbers. 
The definition of the multiple (double in this case) zeta value C(°> ^) is given by 

C(a,b)^^2^—b- ( 8 - 2 °) 
^-^ n a m° 

n>m 

For integer values of a, b this can generically be written in terms of regular Q values and 
products thereof. A useful program on the web for doing this is EZFace 21 . These multiple 
zeta values can then be related to sums of polygamma function via harmonic sums [69] 

£ = (-^WICWCP + 1) - C(«, b + l)]. (8.21) 

0=1 ^ 

We will also encounter more generic sums including ^-functions. Following [5] we can write 
for sums of polygamma's 

f>(Q)V^(Q) = (-l) fe f°° dtC-\Y,)Y(b + 1} L WO ~ C(& + 1) ? (g22) 

where £ _1 (E) is the inverse Laplace transform of E(Q). For the digamma function one 
finds 

f;S(Q)^(Q)= r^-HS) l0g(l ~ e "? +7 

In general when one encounters rational functions Y1q=i ^{Q) 1 P iQ): the recipe that 
proved useful in explicit computations was to split it into two parts E(Q) = Eo(Q) + Ei(Q), 
where 



n>l ^ 



To the term containing Eo one can then apply formula ( |8.21| ) and on the remainder Ei one 
can use the integral relations ( 8.22j ) and ( 3.23] ). 



Let us now apply the above discussion to the integral E^ Q . The integral is most 
easily evaluated by summing residues. The integrand has three different types of poles in 
the upper half-plane 



1 http:/ /oldweb.cecm.sfu.ca/projects/EZFace/ 
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• dynamical poles at u + i(Q ± 1), giving a contribution E^-^q dy n 

• poles at i(Q + 2A;) for k £ Z>o, giving -E^lo fc 

• pole at iQ, resulting in E^ Q ^ 

All of the above poles give rise to residues that contribute to the integral and we will discuss 
them separately. 

Dynamical poles 

Let us first consider the residue at the points v = u + i(Q±l). One straightforwardly finds 

t^nlo = ' \ * — — -2- / 

" =u+i(Q+1) ^ {u + i) J {u 2 + lf(2iQ + u + i) 2 (2iQ + u + iy ' 



8iQ(Q + 1) ( H Q _ iu-l + H iu-l 1 — iu + Hciu-l 



8iQ(Q — 1) ( H n iu+\ + Hiu+i n ~\~ H m+i + Htu+i 
rp „ v nlo = V Q ~ — — * — -2~/ /o 9fi \ 

™Sv=u+i(Q-i)^ ( u _i)J( u 2 + l Y(-2Q + iu + iy(2iQ + u-i) J • [ ' 

Both contributions have to be added and then summed over Q. The obtained expression 
naturally splits in two pieces, one which is purely rational and one that also contains terms 
of the form ip(Q + ±^2»). The latter can be summed by making use of the techniques of 



[5] used to derive equations ( p.21 r p.23|) . This yields the following result for the first piece 



p (4),l i(-l)J(u + i)J 



^( J ) (^) (u 2 + 1) y/ J+2 > (^g) ^ (J+1) (jjg) 



J! 4 (J + 2)! (J + l 



.27) 



The remaining rational part can be more straightforwardly summed. Adding up both 
contributions gives 



(4) -8i J u 
E 



NLOA" 2 J (n 2 + l) 4 (n-i) J 



n 2 + 1 ^J+l) (l^u) ^(J) (1^) ^(J-l) (l^m) 



(J + l)! J! 

(8.28) 



In deriving the above results we used the fact that the Bethe root u satisfies the Bethe 

equations = 1- As it is an intermediate result, the quantity -E^lo dyn * s no * rea ^ 

by itself. It will be canceled in the final energy correction as it should because it solely 
consists of ^-functions. 

Poles at i(Q + 2k) 

The integrand also exhibits a pole at i(Q + 2k) due to the harmonic numbers. The residue 
at this point is easily found to be 

res i(Q+2fc) ^ NLO = 4i Y^°(i(Q + 2k), u). (8.29) 
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In other words, this gives a total contribution of 



4to, fc = 4EE y Q LO (^ + 2fc )^)- 



S.30) 



fc=l Q=l 



It is most convenient to first preform the sum over Q. This yields, again on a solution of 
the Bethe equation, 



^y Q LO (i(Q + 2fc),n) = 
Q=l 

2 w (-i) J u + i-2ik ^ 2 3 ~ 2J i -0( J )(fc+l) 

(t| + i )J (u 2 + 1) S fc J+2 + 2_, " (J) " (n 2 + 1) 3 H+1J, + 



(8.31) 



i4-2J 



U 



{u 2 + If k J [{2k + iu) 2 - 1} 
where 



(u 2 + i)^( J+1 )(fc + i) ^ J \k + i) v (J_1) ( fc + 1 ) 



4m(J + 1)! 



+ 



J! 



22J-n-l 



1 

(« + *) 



n+1 



(«-*) 



u 2 + 1 4ra 



fc J+2 fc J+1 fc J 



m(J-l)! 

^(J-n-l)(fc + l) 
(n 2 + l) 3 r(J-n) 



By using the relation 



^( n )(Jfe + l) = ^ n \k) 



-l) n n\ 
k n+1 ' 



.32) 



.33) 



it is now straightforward to sum first line over k by means of the previously discussed 
techniques. 

The terms from the second line are all proportional to y^k+iu) 2 -!} ' This function 
can easily be split up in the following way 



(2iy 



i 



J-2 



k J [(2k + iu) 2 - 1] {u + i) J (2k + iu) 2 



n=0 



1 



1 



(u + i) n (u-i) n 



(2i) 



n-2 



134) 



The tail allows for a direct expression in terms of multiple zeta values. The first piece has 
to be inverse Laplace transformed 



-i 



1 



1 



(2k + iu) 2 -l) 4 
One can then explicitly preform the integral (|8.22j) to find 



135) 



i 



(2/c + iu) 2 - 1 

2i(-l) fc fc!C(A: + 1) 
it — i 



r(fe + i 



Li fe+1 (e-*)-C(A; + l) 



-l) fc (2i) fc+1 A:! 



e* - 1 
(u + i) k+l 



.36) 



(-l) fc fc!^ 



n=l 



(2i) n ((k - n + 2) 
(it + i) n 
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Pole at iQ 



The residue at v = iQ is the most complicated one. It receives contributions from both the 
rational part and the harmonic number part. Let us expand the harmonic number part 
around this point 



H Q + iv + Hiv-Q + Hg-iv + H Q + iv 



2i 



v — iQ 



n=2 



1 \ n 2(-i) n +V (n) ( 1 ) - V (n) (<2) - 4> {n) {Q + 1) 



(8.37) 
(v-iQ) n . 



For any value of J one can now easily extract the residue at this point and apply the 
discussed techniques. 



8.4.2 Some integrals 

In this section we present a list with integrals that are encountered throughout the main 
text. These integrals are 



1 °° f v — u 

^ = to S J dV Q*-l+l*-u* Y ° Q 

V — 1 



(«), 



_ (u 2 + If 

-htj 



oo „ 



Av 2 



dv 



(Q 2 + v 2 ) 2 
(Q 2 -i 



(Q 2 + v 2 ){\ - Q 2 - v 2 + v 2 ) 



(8.38) 
(8.39) 
(8.40) 



Integral X\ is computed similarly to E-^q. One finds 

A X {J + 1) + B x { J + 1) A x ( J + 2) + B x { J + 2) 



27TlliJ) - (u 2 + m.j+i)\ 

where the coefficients in the above expression are 



( u 2 + i)(j + 2)! 



MJ) = J 

2vri 2 J 



(8.41) 



1 



1 



+ 



(u + iy- 1 (u-iy- 1 
j(j-i) 



(j-i) 



1 — iu 



1 + iu 



+ 



.42) 



2 J - X |_(u + f) J (w-i) J _ 



^-2)(l^<H > ) +( _ 1) ^(J-2)/ r l + ''' 



^(J) _ 2J ^ r(J - rc)r (J - n + \) 



2m 



n=l 



r(J-2n + l) |_(u-i) 2n (u + i) 2 ™_ 



C(2J-2n- 1). 



(8.43) 
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One can then express the remaining integrals 2^3 in terms of the previously computed 
integrals and derivatives thereof 

32£[ Q (J + 2) _ d (u* + l?P£ {J) 
21 ' (J + 2)(J + 3) du J + 3 

^ 4(n 2 + l) 2 7; 1 (J) 2(^ + l) 2 ^ Q (J + l) 

du (J + 2)(J + 3) J + 3 ' 1 j 

1*7) - uX l(J + 1) - + 1) - 4 ^ (J ; J %^ ( ^ )(M ' + 1)1 • (8-45) 
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